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Sunnary 


This  program  was  undertaken  because  of  some  relatively  divergent 
results  in  the  literature  concerning  combustion  noise  from  gas  turbine 
combustors.  Especially  in  regard  to  sound  power  output  the  results  of 
several  workers  could  not  be  reconciled  and  there  was  question  con- 
cerning the  actual  sources  of  noise  which  were  present  in  various  ex- 
periments, The  purposes  here  were  unification  of  past  results,  de- 
velopment of  a rational  foundation  for  the  results  and  noise  source 
identification, 

A theory  was  developed  here  for  the  behavior  of  the  combustor 
interior  acoustics  which  clearly  shows  that  the  interior  pressure 
fluctuations  are  related  to  at  least  three  causal  mechanisms  - direct 
combustion  noise,  entropy  noise  and  vorticity  - nozzle  interaction 
noise.  After  correcting  the  results  of  many  past  experiments  through 
use  of  the  theory  it  is  shown  that  a unified  correlation  of  sound 
power  output  from  combustor  rigs  can  indeed  be  developed,  and,  more- 
over, that  it  has  a rational  theoretical  basis.  It  was  shown  that 
most  past  tests  were  in  configurations  which  precluded  appearance  of 
entropy  and  vorticity  noise.  Theoretically,  however,  it  was  shown  that 
there  is  a reasonable  expectation  that  if  the  level  of  these  two  sources 
could  be  predicted  or  measured  the  scaling  could  be  expected  to  follow 
the  same  laws  as  direct  combustion  noise. 

Experimentally,  the  situation  found  in  engines  was  simulated  - 
the  combustor  terminhted  by  a choked  nozzle.  Here,  vorticity  and  en- 
tropy noise  would  be  expected  to  appear  and,  indeed,  were  found.  By 


use  of  cross  spectral  density  and  multiple  coherence  techniques  It 
was  found  for  the  combustor  tested  that  a)  vortlclty  noise  was  at 
least  as  large  as  combustion  noise,  b)  entropy  noise  Is  minor  Insofar 
as  overall  sound  output  Is  concerned  and  appears  mainly  at  low  fre- 
quency, c)  all  noise  sources  are  coherent  with  each  other  over  various 
frequency  ranges  so  that  an  unambiguous  separation  could  not  be  made, 
d)  due  to  the  choked  nozzle  Impedance  and  transmittance  properties 
the  noise  Is  higher  frequency  2000  Hz)  than  heretofore  believed 
and  e)  some  other  unexplained  noise  source  or  nonlinear  behavior  con- 
taminated the  results  at  low  frequency. 

It  was  additionally  found  that  the  overall  noise  output  was 
higher  than  expected  In  the  configuration  used,  which  employed  a dif- 
fuser terminating  the  choked  nozzle.  It  Is  conjectured  that  the  addi- 
tional noise,  which  was  coherent  with  combustor  Interior  events,  was 
due  to  vortlclty  Interaction  with  the  decelerating  diffuser  flow.  It 
was  also  experimentally  found  that  the  theoretical  expectation  was 
borne  out  that  the  total  noise  scaled  In  accordance  with  the  direct 
combustion  noise  scaling  laws. 

Some  effort  was  expended  In  direct  measurement  of  the  nozzle 
Impedance  properties  through  use  of  a separate  entropy  tube  facility 
and  the  use  of  cooled  film  probes  In  the  combustor.  These  efforts 
were  of  marginal  success.  In  the  end,  theoretical  Impedance  values 
were  used  for  the  source  Identification  portion  of  the  program. 

The  correlation  techniques  of  Pratt  & Whitney  and  General  Electric 
for  combustor  rigs  and  engines  were  reviewed  and  some  suggestions  for 
Incorporating  the  results  of  this  program  were  made. 
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I . Introduction 

During  the  past  decade  there  has  been  rather  intense  research 
on  the  question  of  "core  engine",  "core",  or  "excess"  noise  from 
aircraft  gas  turbine  units.  This  is  noise  which  emanates  from  the 
core  (working  part  of  the  airflow  which  passes  through  the  turbine) 
of  the  engine.  All  noise  due  to  the  fan  (on  a turbofan  engine)  or 

[ 

[ 

^ the  jet  exhaust  is  excluded;  moreover,  turbine  noise,  which  is 

identifiable  by  itself,  is  excluded.  For  several  reasons,  it  has 
been  suspected  for  some  time  that  processes  taking  place  within 
the  combustor  are  responsible  for  core  noise,  but  there  has  been 
■ controversy  over  the  exact  nature  of  the  processes. 

Pratt  AWhitney^^^  claim  correlation  of  combustor  rig  and 
core  noise  data  on  the  basis  of  direct  combustion  noise  as  the 

I dominant  noise  source.  Direct  combustion  noise  is  a dilation  of 

I 

the  flow,  and  consequent  acoustic  wave  production,  due  to  fluctuations 

I in  the  aggregate  heat  release  rate  in  the  combustor.  On  the  other  hand, 

I f 2^ 

► Cumpsty'  ' claims  correlation  of  core  noise  results  based  on  the 

j mechanism  of  entropy  noise.  Entropy  noise  is  noise  produced  by 

I 

I hot  or  cold  spots  traversing  a pressure  gradient  such  as  exists 

’ in  the  nozzle  guide  vanes  preceedlng  the  turbine.  Moreover, 

(3) 

combustor  rig  tests  on  a single  combustor  have  recently  shown 
that  under  conditions  of  such  a pressure  gradient  entropy  noise  tends 

\ 

to  dominate  direct  combustion  noise.  Calculations  in  Ref.  (1), 
however,  based  upon  a particular  theory  of  entropy  noise,  do  not 
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support  the  contention  that  entropy  noise  is  a strong  noise  source 
in  engines . 

These  controversies  are  compounded  by  the  facts  that  recent  com- 
bustor rig  and  core  noise  data,  correlated  without  any  theoretical 

(41 

basis,  for  General  Electric  devices  do  not  agree  with  the 
correlations  of  Kef.  (1).  It  has  furthermore  been  noted  in  Ref.  (3) 
that  there  may  be  another  noise  source  which  has  not  been  considered 
or  measured  as  yet.  This  new  source,  called  vortlcity  noise,  is 
caused  by  mass  flow  fluctuations,  due  to  the  turbulence,  entering 
the  choked  or  nearly  choked  flow  of  the  nozzle  guide  vanes.  Order 
of  magnitude  calculations  suggest  that  this  source  is  at  least  as 
large  as  entropy  noise. 

Consequently,  this  program  was  undertaken  to  see  if  a)  the 
results  of  Refs.  (1)  and  (4)  are,  in  fact,  compatable,  b)  the  three 
sources  of  vortlcity,  entropy  and  direct  combustion  noise  can  be 
separated  and  quantified  and  c)  a unified  core  noise  correlation 
could  be  developed. 


I 
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II.  Theory 


General  Formulation 


The  general  aeroacoustlcs  formalism  to  be  used  for  the  theory 
of  combustor  noise  is  the  linearized  vorticity-acoustic  field 
approach.  This  was  used  in  Ref.  (3)  but  some  of  the  details  were 
omitted  there  and  some  modifications  are  necessary  to  account  for 
vorticity  noise.  The  approach  has  a particular  simplicity  (and  is, 
in  fact,  only  valid)  for  low  frequency  noise  so  that  the  plane  wave 
mode  is  the  only  important  acoustic  mode.  The  approach  is  detailed 
in  Ref.  (5)  for  a constant  property  mean  flow,  but  substantial  modi- 
fications are  necessary  in  the  case  of  combustors  with  variable 
velocity  and  thermodynamic  properties. 

The  approach  consists  of  letting  the  perturbation  velocity, 
and  the  perturbation  pressure,  p^,  be  split  into  a vortical 
part,  associated  with  the  turbulence,  and  a dilational  part, 
associated  with  the  acoustic  wave  motion.  That  is 


V + V 

a 


P + P 
V a 


where  the  velocity  components  have  the  property  that 
Vxv^=V»v^*=0.  These  forms  are  placed  in  the  inviscid 
equations  of  gasdynamlcs  with  chemical  heat  release  and  then  split 
into  a set  of  a)  steady  state  equations,  b)  perturbation  equations 
for  the  acoustic  components  and  c)  perturbation  equations  for  the 


4 


vortical  components.  This  Is  done  In  Ref.  (3).  However,  It  Is 
found  that  such  a spilt  Is  only  rigorously  possible  If  the  condi- 
tion 


Vp^x^p*0 

A 


Is  met.  That  Is,  the  gradient  of  the  acoustic  pressure  must  be 

parallel  to  the  gradient  of  mean  density.  The  only  reasonable  flow 

and  acoustic  situation  which  will  satisfy  this  requirement  Is  a 

nearly  one-dlmenslonal  mean  flow  upon  which  Is  superimposed  plane 

wave  acoustic  motion.  This  will  be  assumed. 

The  limit  of  low  Mach  number  will  also  be  assumed.  In  which 

case  the  differential  equations  for  the  Fourier  transform  of  the 

(3) 

acoustic  pressure  becomes 


+ k 


^U) 


f 


(U 


(V  - 1) 


(1) 


Also  assumed  In  £q.  (1)  Is  that  the  speed  of  sound  Is  roughly  con- 
stant through  the  combustor  (an  average  value  has  been  used) . The 
"a"  subscript  has  been  dropped  for  clarity  In  Eq.  (1)  but  It  Is 
understood  that  the  acoustic  pressure  obeys  Eq.  (1).  Also,  to  be 
consistent  with  the  splitting  requirement  the  transverse  derivative 
In  Eq.  (1)  must  be  weak  compared  with  the  longitudinal  derivative. 

The  neglect  of  terms  proportional  to  the  Mach  number  In  Eq.  (1) 
la  permlssable  because  there  are  no  essential  physics  omitted. 

That  Is,  the  effect  of  the  flow  terms  Is  minor  and  Is  primarily 
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( 

1 
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felt  In  the  wave  propagation  speed  so  that  resonant  frequencies  are 
only  subtly  shifted  from  their  no-flow  values.  This  Is  also  the 
primary  effect  of  variable  speed  of  sound.  The  same  is  not  true  when 
one  considers  the  boundary  conditions  to  Eq.  (1),  where  flow  terms 
are  essential  to  a proper  description  of  acoustic  damping  and 
noise  generation,  as  will  be  seen.  Consequently,  more  care  Is 
necessary  In  the  boundary  conditions  formulation.  It  has  further- 
more been  assumed  In  Eq.  (1)  that  llquld-gas  momentum  (drag) 
Interactions  are  unimportant.  This  has  justification  In  numerous 
studies  of  combustion  instability, where  the  drag  effect  is  of 
the  same  order  of  magnitude  as  the  flow  terms. 

Cross-section  averages  of  terms  will  be  denoted  by  a hat  super- 
script. That  is,  for  example, 

p = f dA  p (r)  = p (x) 

f - 7 f dA  f (s) 

(1)  A «)  u)  “• 

These  are  the  effective  plane  wave  values  of  quantities.  Taking 
the  cross-section  average  of  Eq.  (1)  and  assuming  a constant  area 
combustor 


P, 


dx 


^ + p 


uu 


H 


^p 


(U 


n dC 


- f 


U) 


(2) 


f 

I 

I 

:i 

ii 
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where  the  divergence  theorem  has  been  used. 
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P which  will  later  be  experimentally  determined,  but  contains 

w 

the  wall  damping  effect.  In  Eq.  (4),  therefore,  the  wall  boundary 
condition  is  already  Incorporated. 

At  the  combustor  head  end  it  is  assumed  the  wall  (even  if 
there  are  air  entrance  slots)  is  acoustically  hard  so  that 
u^(0,y,z)  » u^(0)  = 0.  Equation  (3)  then  yields 


) ■» 

x=0 


At  the  nozzle  end  an  admittance  condition  is  involved  for 
the  plane  wave  quantities  that 


+ b Z = 0 
m (u 


In  several  limits  the  admittance  coefficients,  a^  and  b^  may  be 

calculated.  For  example,  if  the  nozzle  is  choked  and  the  frequency 

is  low^^^  a - (y-1)  M /2v  and  b - -M-/2  ; if  the  nozzle  is 

a pipe,  open  to  free  space,  and  the  Mach  number  is  low,  the  Levine- 

Schwinger  problem is  the  appropriate  limit,  whereby  as  frequency 

tends  toward  zero  a^  “•  “ and  b^  = 0.  In  general,  however,  a^ 

and  b must  be  measured. 

U) 

Taking  the  axial  component  of  Eq.  (3)  and  averaging  over  the 
cross  section, 
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L 


written  in  the  form  of  the  second  law  of  thermodynamics,  demands 
that  the  entropy  fluctuations  be  convected  with  the  fluid.  Hence, 


d S 


M -7-^  + i k E =0 
e dx  u) 


at  the  exit  plane  and  Eq.  (10)  simplifies  to 


dx 


uja 


M + 


^ C 
Ya  +7 
U)  A 


M P -) 
e w 1 

ik  J 


u 

r u)  -1 

ikYp  — + b Z (11) 

L u)  (uj  ^ ' 

The  second  term  in  brackets  is  generally  of  the  order  of  and 

would  ordinarily  be  neglected  if  the  procedure  used  to  derive 
Eq.  (1)  were  used.  However,  this  term  provides  damping  and  prevents 
unbounded  resonances.  It  is  important  to  retain  mean  flow-dependent 
terms  when  they  add  danqping,  whereas  in  Eq.  (4)  the  flow  terms  neglected 
only  shift  resonant  frequencies  a slight  amount. 

The  right  hand  side  of  Eq.  (11)  contains  two  source  terms  - 
vorticity  and  entropy  at  the  combustor  termination.  The  combustion 
noise  source  term  occurs  in  the  differential  Eq.  (4).  Damping  is 
provided  by  the  wall  impedance  in  Eqs.  (4)  and  (11),  by  the  nozzle 
through  a^  in  Eq.  (11)  and  by  convection  of  acoustic  energy  by 
the  mean  flow  in  Eq.  (11).  Hence,  all  of  the  basic  physics  are 


10 


Included,  although  some  approximations  have  been  Introduced. 

Solution  and  Orders  of  Magnitude 

Now  dropping  the  "hat"  superscript  and  "a"  subscript;  the 
solution  to  Eq.  (4)  subject  to  Eqs.  (5)  and  (11)  may  readily  be 
obtained.  It  is 


P, 


= A H.  „ + V H„+a  H_ 

U)  AP  U)  vp  (1)  Op 


H 


AP 


H 


vp 


^ i k V iiT/Tg  ^ 


H, 


ap 


cos  K (x-t)  + -jf  sin  K(t-x) 


cos  K X 


U)  V, 

— cos  K X 
M 

e 


0 = 0 cos  Kt  - K sin  Kt- 
e 


P M 3 

a “-ik(M  +'Ya  +—  “ ) 

e e (uAik' 


(12) 


The  three  source  terms,  A , V and  a are  defined  in  the 

U)  (1)  U) 

Nomenclature  and  are  basically,  the  Fourier  transforms  of  the 
fractional  heat  release  rate,  exit  plane  vortlcity  and  exit  plane 
entropy  fluctuation.  From  the  first  of  Eqs.  (12)  the  pressure 
transform  is  linearly  related  to  the  three  causal  processes. 


It  is  desirable  to  gain  an  idea  of  some  orders  of  magnitude 
from  Eqs.  (12).  Assume,  for  the  moment  that  =0  so  that  only 

combustion  noise  is  present.  Furthermore,  assume  that  a combustor  is 
terminated  by  an  impedance  matching  device  so  that  only  downstream 
acoustic  waves  are  propagating  (no  reflections).  In  such  a case 
a » 1/y  and  if  M < < 1 then  a - i k.  Moreover,  for  low 
k « K.  With  these  approximations,  at  x = t. 


P (^) 


Yi  M AT/T 
e e 

- i cos  k-t  + sin  kt 


A 


U) 


Constructing  a spectral  density  from  the  Fourier  transform 


(9) 


2 n 


P,..P 


U)  U) 


(V  M 
P 


AT/Tg)' 


The  acoustic  power  is  given  by 


P * A 


r s 

J 

-oo  p c 
e e 


^e  ■ 2 r“ 

duj  =»  — (Y  M AT/T  ) S. 

- - e e J A 

DC  -<» 

e 


dtt) 


A c ^ p (M  AT/T  )^  < A^  > 
e e e e e 


It  is  to  be  noted  that  a reference  inlet  velocity  may  be  defined 

by  P.  u - = p u and  p./p  = T /f . . Consequently,  the 
■'  i ref  e e ^i  ^e  e i 

grouping  on  the  right  above  may  be  written  in  terms  of  more  con-* 
venient  design  variables  as 


12 


? - MT^2  Y 

^ P ['ZT'i  1 — ^ 

ref  \T^/  R(T^  + AT) 


P = < A > A 

e 


(13) 


Now,  a correlation  formula  is  to  be  presented  later  in  this 
report,  which  accurately  predicts  combustion  noise  experiments  which 
closely  meet  the  requirements  of  the  assumptions  above,  has  been 
applied  to  a particular  run  from  Ref.  (4),  The  numbers  are 


ref 

276  ft/sec 

A = 225  in^ 
e 

p - 65.5  psia 

AT  = 

910  °R 

T.  = 1.635  ‘'r 

1 

P 152.4  w 

2 -4 

When  introduced  into  the  formula  above  < A > = 8.9  x 10  results. 
Consequently  the  rms  value  of  A is  roughly  C.C3.  Another  way  of 
saying  this  is  that  the  fluctuation  level  in  tne  total  heat  release 
ra./a  is  about  3%  of  the  mean  heat  release  rate. 

Now  assume  that  this  heat  release  rate  fluctuation  remains  at  the 
same  level  in  an  installed  conf igiiration  where  the  comb'.istor  is 

terminated  by  a choked  or  nearly  choked  nozzle.  In  such  a case 

M 

~ ® ( 7 ) 

~ for  short  nozzles.  Con.sequent  ly  in  looking  at 

H and  H.  in  Eqs.  (12),  all  three  termt'-  in  the  equation  for 
a p Ap  ■ 

p^/p  scale  almost  exactly  the  same  and  in  fact  are  of  the  same 

2 2 2 

order  of  magnitude  if  < A > , < V > and  < cr  > are  of  the 
same  order  of  magnitude.  This  is  oxpecteA  to  be  so,  since  relative 
intensities  of  turbulence  are  high  in  a combus  tor * and  tempera- 
ture fluctuation  levels  are  also  high,^^*  Of  co^irse  a and  V 
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are  cross-section  averages  which  will  depress  them  in  magnitude 
compared  with  their  point  magnitudes.  But  it  is  expected  that 
these  quantities  will  be  of  the  order  of  1 to  3%.  The  conclusion 
is  that  for  the  interior  pressure  fluctuation,  it  may  be  expected 
to  be  strongly  influenced  by  all  three  mechanisms . 

This  does  not  say  anything  about  the  noise  that  gets  out, 
since  a)  there  are  transmission  properties  of  the  combustor 
termination  to  worry  about  and  b)  the  nozzle  generates  additional 
downstream  noise  by  the  entropy  and  vorticity  mechanisms. 

Using  Equation  (7)  one  can  calculate  the  acoustic  velocity. 
Upon  so  doing  one  discovers  the  highly  important  fact  that  the 
velocity  due  to  the  vorticity  source  is  90°  out  of  phase  with  the 
pressure.  Consequently,  there  is  no  acoustic  termination 
downstream  of  the  combustor  exit  plane  to  interact  with  V to 
produce  noise  no  sound  power  output  can  result  from  V.  Moreover, 


b = 0 in  cases  where  no  pressure  gradient  terminates  the  com- 


U) 


bustor  flow  so  that  these  are  experimental  situations  where  com- 
bustion noise  is  the  only  noise  source. 

Returning  to  Equation  (13),  this  equation  gives  an  idea  of 
scaling  laws  to  be  expected  in  the  absence  of  entropy  and  vorticity 

noise.  The  only  major  problem  is  in  the  estimation  of  the  scaling 
2 


of  < A~  > . For  order  of  magnitude  purposes  the  cos  K x term  in 
the  definition  of  A can  be  omitted  and 


< A^  > 


(’ 

[ q'  dv) 

/ 1 

CM 

1 

1 Q dV^ 
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Estimation  of  the  square  of  stochastic  variables  like  Q'  when 
integrated  over  space  yields 


< A > « V V 


cor 


< > 

av 

(JT^ 


V 

cor 

V 


Mathews assumed  that  the  correlation  volume  was  the  entire 
volume  of  reactants  produced  around  a single  fuel  nozzle.  While 
it  may  be  plausible  that  there  is  a certain  degree  of  correlation 

within  the  volume  served  by  a single  fuel  nozzle,  since  there  are 

I 

; many  uncorrelated  turbulent  eddies  within  such  a volume,  there 

i cannot  be  perfect  correlation.  Moreover,  the  turbulence  size 

I 

[ scale  is  probably  most  closely  linked  to  a transverse  combustor 

dimension,  as  in  pipe  flow.  Consequently,  ‘ 

seems  reasonable  to  assume  with  0 < a < 1.  Then  if  it  is  argued 
that  the  mean  fluctuation  level  of  is  proportional  to  the 

mean  heat  release  rate,  an  assumption  which  frequently  is  borne 
out  as  true  in  turbulence. 


results.  When  this  is  placed  back  into  Eq.  (13) 


„ _ -a  ..  1/2 . 2 - /AT^‘  1 

? « (A  /</)  A u , p [—/  = : — 

f e ' e ref  ^ + AT 


(14) 
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results.  This  suggests  the  variables  against  which  to  test  and 
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correlate  the  data.  Moreover,  it  suggests  that  a power  law  type 
correlation  fit  may  be  adequate. 


Inversion  for  the  Source  Quantities 

The  two  types  of  x-wise  variation  in  the  H's  of  Eqs.  (12) 
suggest  that,  if  a third  independent  measurement  is  made,  two 
different  pressure  measurements  will  allow  an  inversion  to  the 
source  quantities.  In  fact,  taking  positions  x and  x = t and 
denoting  the  pressure  transforms  there  as  p^  and  p^,  there 
results 

i K T (p  cos  K - p cos  K x) 

A = 2 (15) 

V k AT  Mg  sin  ■ x) 


This  highly  important  result  states  that  the  combustion  rate 
fluctuation  may  be  obtained  by  a two  point  pressure  measurement. 
and  Is  independent  of  the  admittance  relations . 

The  third  independent  measurement  which  will  be  made  is  the 
cross  section  average  of  the  temperature  fluctuation  at  the  exit 
plane.  This  is 


T = 
U) 


_1. 

A 

e 


P 


(16) 


It  should  be  noted  that  / and  a are  related  by 

L ta 


I 

ui 


cu 

(AT/T  ) a • Then  Eqs.  (12)  may  also  be  inverted  for  V . 

6 U)  ^ 


The  result  will  not  be  stated  here,  but  V does  depend  upon 
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and  b^,  but  it  is  the  only  one  of  the  three  source  terms 
that  does,  using  this  measurement  scheme. 


It  is  to  be  noted  in  Eq.  (15)  that  the  inversion  process 
becomes  singular  (0/0)  as  u),  k and  K 0.  This  is  the  con- 
sequence of  the  fact  in  Eqs.  (12)  that  there  is  no  wave  structure  as 
frequency  tends  towards  i:ero;  only  a bulk  mode  oscillation  results 
and  = p^.  This  fact  will  be  used  to  advantage  later  in  transducer 
calibration. 

A contaminant,  not  accounted  for  in  the  above  theory,  is 

hydrodynamic  noise.  It  is  non  propagational  pseudosound  which 

sometimes  makes  its  appearance  in  the  combustor  to  be  used  at  low 
(3) 

frequency  . It  can  be  detected,  so  that  if  it  appears  the 
frequency  ranges  where  it  appears  will  be  excluded  from  the  data 
reduction  process.  Otherwise,  it  is  assumed  that  the  theory  is 
valid  and  A,  V and  a may  be  extracted  from  the  internal  measure- 
ments . 


Comparison  between  Experiments  on  Different  Rigs 

Often  in  the  literature  "apples  and  oranges"  have  been 

compared  when  comparing  combustor  noise  results.  The  solution  of 

Eqs.  (12)  clearly  shows  that  the  combustor  termination  (throu^  a ) 

e 

plays  a strong  role  in  the  interior  pressure  fluctuation.  Conse- 
quently, it  may  be  expected  to  play  a strong  role  in  the  noise 
output  also.  Different  terminations  have  been  used  on  different 
experiments.  In  order  to  make  valid  comparisons  between  experiments. 


the  results  must  be  corrected  to  equivalent  terminations. 

Three  cases  will  be  considered.  These  are  a)  a combustor 
terminated  by  an  open  unflanged  end,  open  to  the  atmosphere,  b) 
a termination  by  a short,  choked  nozzle  and  c)  a termination  by 
either  a perfect  impedance  matching  device,  or,  what  is  equivalent, 
an  infinite  extension  to  the  pipe  at  the  exit  plane.  Where  necessary 
it  is  assumed  that  the  combustor  is  a right  circular  cylinder. 

Direct  combustion  noise  is  the  only  noise  source  under  consideration 
and  it  is  assumed  that  is  unaffected  by  the  termination,  for 

the  same  operating  conditions  (there  is  no  reason  to  suspect  other- 
wise). The  exit  Mach  number  is  assumed  low  and  ka  < < 1 so  that 

/o\ 

in  case  a)  the  Levine-Schwinger  result  may  be  applied  that'' 


ou  'Y  k-a  (0.6  i + 0.25  k a) 


In  case  b)  the  nozzle  is  assumed  short  so  that 


a = - ^ M 

U)  2y  e 


In  case  c)  the  Mach  number  is  assumed  low  so  that 


1 

a = - — 

u)  V 

For  purposes  of  calculation  is  assumed  small  so  that  K “*  k. 

Where  necessary  (in  case  a))  ^,/a  = 


7.52  has  been  taken,  correspond- 
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When  only  combustion  noise  Is  under  consideration  the  exit 

plane  pressure,  along  with  the  exit  plane  Impedance,  are  sufficient 

to  calculate  the  power  output.  Using  Eqs.  (12),  the  admittances 

above,  and  calculating  the  power  per  unit  frequency  Interval,  the 

results  may  be  expressed  In  the  form  P /P  and 

(u  case  a u)  case  c 

P . /P  The  results  for  M =0.3  are  shown  In  Figure  1. 

(U  case  b u)  case  c.  e ° 

It  Is  seen  that  the  results  are  quite  frequency  sensitive,  but  that 
either  of  the  reflecting  terminations  causes  a substantial  drop  In 
radiated  power,  compared  with  a non-reflecting  termination.  For 
the  open  end  (case  a)  the  quarter  and  3/4  wave  resonances  are 
clearly  seen,  whereas  for  the  short  choked  nozzle  case  the  half-wave 
resonance  Is  present.  This  result  shows  clearly  that.  If  results 
from  case  c and  case  a,  for  example,  are  compared  experimentally,  a 
correction  of  the  order  of  10  dB  In  radiated  power  must  be  applied  to 
the  results.  Moreover,  the  correction  must  carefully  consider  the 
frequency  content  of  the  noise.  Such  a correction  will  be  used  in 
the  next  section. 

Exterior  Noise 

The  transform  of  the  exterior  noise  may  always  be  written  as 


H. , A + H , a + H,,,  V + Z 
A1  (u  cl  U)  vl  U)  a 


where  Z Is  the  deviation  from  the  linear  causality  of  the 


Figure  1,  Theoretical  comparison  of  the  sound  radiation 
per  unit  frequency  interval  for  two  different 
exit  plane  terminations.  Only  direct  combustion 
noise  is  considered. 
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first  three  terms.  It  is  known  that  even  if  Z is  substantial 

(U 


that  the  best  linear  fit  in  a least  squares  sense  is  obtained  by 
assuming  Z^  = 0 and  measuring  the  H's  in  Eq.  (17).  This  is  done 
by  constructing  cross  spectral  densities  in  the  form 

SlA  = »AlSAA^«al|A  + «Vl  ^VA 

= »A1  ^Aa  »al  ^ ^1  \ 

» 1 "aV  . (18) 

measuring  the  cross  spectra  and  solving  Eqs.  (18)  for  the  transfer 
functions.  In  practice,  the  cross  spectra  Sj^q, 

and  are  measured  and  the  cross  spectra  of  the  right 

hand  sides  of  Eqs.  (18)  are  calculated  by  inverting  Eqs.  (12)- 
that  is  ,through  use  of  Eqs.  (15)  and  (16)  and  the  associated 
equation  for  V^. 

The  goodness  of  the  fit,  or  stated  otherwise,  the  magnitude 

of  the  contaminant  Z^  ,1s  determined  by  the  multiple  coherence 
(9'i 

function'  '' 


''ll  ■ (“ai  Su  * "vi  hv  * =1,  )'®11 
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Therefore  multiplication  of  the  measured  spectrum,  by 

gives  that  portion  of  the  spectrum  which  is  linearly,  causally 
related  to  the  measured  sources.  It  is  called  the  multiple  coherent 
spectrum. 

If  the  off  diagonal  elements  of  Eqs.  (18)  are  much  smaller 
in  magnitude  than  the  diagonal  elements,  the  three  sources  are 
uncorrelated.  This  would  be  the  ideal  situation  since  a definite 
splltof  the  noise  into  three  components  would  be  possible.  It  will 
be  found,  however,  that  this  did  not  occur  in  certain  frequency 
banas.  A measure  of  correlation  between  the  sources  is  the 
ordinary  coherence  function  given  by 


iil 


s,,  s.. 

ii  JJ 


0 i Y s 1 


This  function  will  also  be  used  in  the  following  sections. 
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III.  Correlation  of  Past  Combustor  Rig  Tests 

In  References  (1),  (4),  (10)  and  (11)-(12)  there  is  a 
wealth  of  data  regarding  combustion  noise  from  combustor  rig  tests. 
Virtually  all  relevant  parameters  have  been  varied, tests  at  simulated 
engine  conditions  have  been  conducted,  and  different  combustor  types 
were  used.  Comparison  of  correlations  which  have  been  produced  for 
the  data,  however,  shows  that  the  correlation  of  one  manufacturer 
does  not  recover  that  of  another  manufacturer.  It  is  the  purpose 
here  to  see  if  a unified  correlation  can  be  developed  which  will 
adequately  cover  the  data  ranges  of  the  cited  references. 

The  work  here  will  concentrate  on  sound  power  only.  The 

frequency  content  of  the  cited  references  poses  one  of  the  great 

(13) 

unsolved  problems  in  combustion  noise.  In  Reference  (1)  a 

theory  is  presented  for  the  frequency  content  which  appears  to 
recover  the  correct  trends  for  the  combustors  of  that  work,  but  it 
works  for  no  other  manufacturers Moreover , the  results  of  Ref.  (10) 
are  not  explained  by  that  theory.  About  all  that  can  really  be 
said  is  that  combustion  noise  has  been  observed  in  the  past  to  be 
basically  low  frequency  in  nature  (^  1000  Hz).  Of  course,  it  must 
be  borne  in  mind  that  the  frequency  content  depends  upon  the 
acoustics  of  the  containing  hardware.  Since  the  combustor  rigs 
have  been  different,  it  is  perhaps  not  surprising  that  there  is 
confusion  concerning  the  observed  frequency  content.  In  any  event, 
even  in  this  program  there  will  be  found  some  new  twists  in 
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frequency  content  so  that  the  issue  will  be  avoided  here.  Only  sound 
power  will  be  considered. 

In  the  rigs  of  Refs.  (1)  and  (10)  the  six  different  combustors 
were  exhausting  directly  (with  some  diffusion  before  exit)  to  the 
atmosphere.  The  acoustic  configuration  was  essentially  that  of  an 
open-ended  organ  pipe  and  resonances  are  clearly  seen  in  the  data  at 

* 

the  1/4,  3/4,  5/4,  etc.  wave  resonance  points.  On  the  other  hand,  the 
work  of  Refs.  (4)  and  (11) -(12)  was  carried  out  in  two  different  rigs 
which  were  assumed  reflection  free  at  the  combustor  termination.  The 
data  reduction  procedure  assumed  that  waves  emanating  from  the  com- 
bustor were  not  reflected  from  downstream  impedance  mismatches.  In 
viewing  these  rigs,  it  is  the  judgement  here  that  such  an  assumption 
was  probably  justified  at  least  to  a reasonable  approximation.  Con- 
sequently, the  data  of  the  various  rigs  cannot  be  directly  compared 
unless  corrections  are  applied  to  bring  the  data  to  the  same  basis, 
Hie  basis  which  is  chosen  here  is  to  correct  all  acoustic  power 
measurements  to  the  "equivalent  reflection  free  power",  Thus, 

the  work  of  Ref,  (4)  and  (11)-(12)  will  be  accepted  as  island  correc- 
tions will  be  applied  to  the  data  of  Refs.  (1)  and  (10).  'Hiis  cor- 
rection will  be  simply  done.  That  is,  usually  in  Ref,  (1)  the  com- 
bustion noise  spectrum  peaked  at  a frequency  slightly  higher  than 
the  3/4  wave  resonance  point.  Reference  to  Fig,  1 shows  that  a cor- 
< rection  of  + 7 dB  is  a reasonable  correction  to  add  to  the  data  of 


Ref,  (1),  In  Ref,  (10)  most  of  the  noise  was  centered  about  the  1/4 
wave  resonance  point,  and  reference  to  Fig,  1 shows  that  a correction 


of  + 1 dB  is  reasonable.  For  high  accuracy  a frequency  dependent 
correction  should  be  applied  and  a summation  over  frequency  bands 
employed.  However,  it  will  be  seen  that  the  crude  approximation 
above  yields  satisfactory  results. 

The  next  thing  to  notice  is  that  no  entropy  or  vorticity  noise 
is  expected  in  the  cited  rig  tests,  as  discussed  in  Section  II,  Con- 
sequently, there  is  high  confidence  that  a form  of  the  correlation 
law  suggested  by  Eq,  (14)  will  satisfactorily  correlate  the  data. 

The  form  chosen  is 

P = a,  p^2  u /3  a4  ^aj  ^ ay 

eq  1 ref  i f e ^ e 

(21) 

This  form  is  close  to  the  form  of  Eq,  (14)  and  is  suitable  for 
linear  multiple  regression  analyals  after  taking  the  logarithm  of 
Eq.  (21),  202  data  points  were  available  from  the  cited  references, 
with  15  different  combustors.  The  regression  analysis  yielded  the 
following  results,  where  the  theoretically  expected  values  are  also 
shown: 

Table  1 

Regression  Analysis  Results  and  Comparisons  with  Theory 


*1 

*2 

^4 

®5 

^6 

Experiment 

0.91 

1.9 

3.4 

-2.5 

1.3 

-0.78 

1.0 

1.0 

Theory 

N/A 

1 

2 

-2  to 

-3 

2 

0 to 
-1 

1 

1 

Note: 

The 

units  are 

power 

in  watts. 

pressure 

in  psia. 

velocity  in  ft/sec,  inlet  temperature  in  °R  and 
exit  area  in  square  inches.  All  other  quantities 
are  dimensionless 
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The  standard  deviation  of  the  correlation  is  + 3 dB  and  a comparison 
of  the  regression  formula  with  the  data  is  shown  in  Figo  2.  Ihe  in- 
dividual General  Electric  combustors  are  pointed  out  in  this  figure 
because  it  was  with  these  combustors  that  the  largest  scatter  occurred. 
The  correlation  is  considered  highly  interesting  for  several  reasons. 
These  are  a)  there  is  a rational  theoretical  basis  which  generally 
holds  up,  b)  the  single  correlation  formula  holds  over  an  extreme 
range  of  acoustic  power  (10^)  and  c)  the  Ref,  (1)  conclusion  of  a 
fuel  nozzle  effect  holds  up  when  all  combustors  are  considered. 

A close  look  at  the  General  Electric  designs  revealed  no  design 
parameter  which,  if  included  in  the  correlation,  would  reduce  the 
scatter.  Consequently,  the  + 3 dB  is  accepted  as  mostly  due  to  experi- 
mental error  in  the  variety  of  rigs  employed  and  in  data  analysis 
methods. 

It  is  to  be  noted  that  if  the  results  of  Ref.  (1)  are  not  changed 
by  the  7 dB  correction,  the  a^  exponent  comes  out  as  2,0  rather  than 
1,0,  This  is  not  acceptable  from  the  theoretical  viewpoint  and  is 
viewed  as  confirmation  of  the  correction  procedure. 

The  primary  deviation  between  theory  and  experiment  lies  in  the 
exponents.  There  is  an  interesting  interpretation  for  this. 
The  deviation  from  theory  is  by  a factor  of 
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That  is,  the  deviation  looks  closely  related  to  the  dynamic  head,  a 
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Figure  2.  Comparison  of  the  predicted  power  from  the  re- 
gression equation  with  the  experimental  sound 
power. 
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a term  which  often  enters  in  aerodynamic  noise  computationso  Moreover, 
aerodynamic  noise  usually  scales  with  a high  exponent  on  velocity. 

It  is  conjectured  here  that  since  no  experiment  can  be  entirely  free 
of  aerodynamic  noise  contamination  of  the  combustion  noise  results 
the  high  and  a^  exponents  may  be  false.  It  is  possible  that  the 
theory  is  more  accurate  as  to  combustion  noise  scaling  laws  than  the 
experiments. 

It  is  now  emphasized  that  this  is  a combustor  rig  correlation 
for  equivalent  power  in  a reflection  free  situation  and  in  the  absence 
of  entropy  and  vorticity  noise.  In  an  engine  installation  where  the  com- 
bustor is  terminated  by  nozzle  guide  vanes,  the  turbine  assembly  and 
exhaust  duct, several  things  change  tdiich  will  change  the  acoustic 
power  output.  First,  the  nozzle  guide  vanes  usually  operate  at  a 
choked  or  nearly  choked  condition.  This  necessitates  corrections  by 
a form  of  Fig,  1 and  also  brings  in  the  possibility  of  entropy  and 
vorticity  noise.  Second,  more  noise  due  to  vorticity  and  entropy  can 
be  created  in  the  turbine,  and  this  is  combustor-related.  Finally, 
there  are  attenuation  processes  in  the  turbine  and  nozzle  which  must 
be  accounted  for.  These  issues  will  be  addressed  in  Section  VIII. 
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IV,  Combustor  Apparatus  and  Experimental  Procedures 
Combustor  Details 

The  combustor,  used  in  the  present  investigation,  is  of  can- 
type,  taken  from  a Boeing  502-7D  gas  turbine  engine  unit.  The  origi- 
nal air  inlet  has  been  modified  to  suit  the  experimental  set-up.  Fig- 
ure .3  shows  a cross  sectional  view  of  the  combustor.  Liquid  fuel  is 
sprayed  into  the  burner  liner  through  a fuel  nozzle.  Air  for  com- 
bustion enters  through  the  slots  in  the  head  end  and  walls  of  the 
liner.  The  by-passed  air  leaves  the  combustor  through  the  annular 
space  between  the  liner  and  the  burner  shell.  The  geometry  and  the 
arrangement  of  the  slots  is  shown  in  Fig,  3,  A spark  plug  is  used  to 

ignite  the  fuel/air  mixture.  Stable  operation  limits  for  this  com-  ' 

bustor  were  determined  in  a previous  program  (Reference  (10)),  The  . 

i 

by-pass  air  flow  rate  is  measured  to  be  26%  of  total  air  supplied, 

Pitot  pressure  distribution  and  mean  temperature  distribution  of  the  ^ 

hot  gases  at  the  burner  exit  plane  are  reported  in  Ref.  (3). 

1 

Flow  Systems 


The  schematic  lay-out  of  the  flow  system  is  shown  in  Fig,  4,  ! 

. ♦ ' 

Air  for  the  combustor  is  supplied  from  an  861  KPa  air  reservior, 
regulated  by  a valve  and  metered  by  an  orifice  meter,  A muffler  is 

I 

provided  in  the  line,  upstream  of  the  combustor,  to  reduce  flow  I 


Figure  3,  Combustor  Details. 


, Flow  System  Schematic. 


tor  are  all  mounted  on  a test  stand  outside  the  laboratory  at  an 
open  site  while  the  flow  control  valves,  the  flow  meter  displays 
and  the  data  acquistlon  systems  are  situated  inside  the  laboratory 
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Measurement  Techniques 

The  main  objective  of  this  experimental  program  is  to  estimate 
the  strengths  of  various  noise  sources  located  inside  the  combustor 
(viz,  combustion,  entropy  6e  vorticity)  and  their  relative  contribu- 
tions to  the  exterior  radiated  noise.  This  requires  measurements  of 
interior  pressure,  temperature  and  velocity  fluctuations  as  well  as 
exterior  noise  characteristics.  Moreover,  it  was  desired  to  simulate 
an  engine  condition  of  choked  nozzle  guide  vanes.  Consequently,  the 
combustor  under  consideration  was  fitted  with  three  terminating  con- 
verging-diverging nozzles.  These  provided  contraction  ratios  of  14.6, 
12,1  and  9,1,  The  diverging  section  is  essential  to  keep  jet  noise 
levels  very  low,  thereby  eliminating  this  contaminant  from  the  ex- 
terior measurements.  After  the  first  test  run,  it  was  discovered  that 
a flow  separation  was  occuring  upstream  of  the  diffuser  exit.  This  was 
due  to  a faulty  assumption  of  considering  the  diffuser  as  two-dimen- 
sional rather  than  three-dimensional.  This  resulted  in  too  large  an 
expansion  angle.  A new  diffuser  was  built  with  very  low  expansion 
angle  (half  angle  = 2°)  to  remedy  the  separation  problem.  For  pri- 
mary data  reduction,  four  tests  were  carried  out,  all  at  choked  con- 
ditions. Three  of  them  \rere  conducted  with  the  above  three  nozzle- 


diffuser  assemblies.  Another  one  was  at  the  middle  contraction  ratio. 
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with  only  the  nozzle  without  diffuser  in  place.  The  operating  con- 
ditions for  these  four  tests  are  listed  in  Table  2,  Apart  from  these 
main  tests,  for  the  middle  contraction  ratio,  airflow  and  fuel/air 
ratio  traverses  were  run  a)  at  constant  airflow  with  varying  fuel/air 
ratios  and  b)  at  constant  fuel/air  ratio  with  varying  airflows. 

Pressure  Fluctuation  Measurements 

The  combustor  interior  pressure  fluctuations  are  measured  by 

Fhotocon  transducers,  mounted  at  locations  slightly  displaced  from 

the  flow  field  through  an  infinite  tube  assembly.  This  arrangement  is 

necessary  to  eliminate  contamination  by  hydrodynamic  (turbulence) 

(3) 

noise.  The  photocon  transducer  can  extract  signals  above  130  dB 
-5  2 

referred  to  2 x 10  N/m  . The  interior  pressure  measurements  are 
made  at  two  points-  one  downstream  of  the  primary  combustion  zone 
but  upstream  of  the  nozzle  and  the  other  at  the  nozzle  entrance  plane 
as  shown  in  Figure  5. 

The  exterior  radiated  noise  is  measured  by  three  near  field 
BiWel  £c  Kjaer  type  413A,  12  mm  condenser  microphones  at  45°,  90°  and 
135°  to  Jet  axis.  All  the  three  microphones  are  situated  at  a radius 
25.5  cm  from  the  combustor  exit  plane.  The  reason  for  choosing  near 
rather  than  far  field  is  that  the  latter  measurements  are  contaminated 
by  reflections  from  ground  and  nearby  constructions.  Ihe  relative  po- 
sitions of  the  interior  and  exterior  pressure  transducers  are  shown 
in  Figure  5.  All  these  transducers  are  calibrated  by  a VAilttaker  type 
PC-125  acoustic  calibrator.  The  near  field  pressure  measurements  give 
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an  accurate  measurement  of  acoustic  power  because  of  the  basically 

(3) 

monopole  radiation  pattern. 

Temperature  Fluctuation  Measurements 

The  entropy  noise  study  requires  the  details  of  combustion 
generated  hot  spots;  hence,  the  measurement  of  temperature  fluctua- 
tions of  hot  gases  at  the  burner  exit  (nozzle  entrance)  plane  have 
been  carried  out  using  fast  response  Chromel-Alumel  thermocouples. 

In  the  present  program,  a thermocouple  junction  constructed  from  a 
wire  diameter  of  0.025  mm  is  selected  as  a compromise  between  the 
conflicting  requirements  of  longevity  and  response  times. 

When  a thermocouple  is  used  as  the  primary  element  for  rapidly 
fluctuating  temperature  measurements,  the  thermal  lag  of  the  element, 
caused  by  the  slow  rate  of  heat  transfer  between  the  element  and  the 
surrounding  gas,  limits  the  ability  of  the  element  to  follow  such 
fluctuations  truly.  This  thermal  lag  is  quantified  by  a parameter 
called  the  "time  constant"  of  the  thermocouple.  The  value  of  the  time 
constant  depends  mostly  on  the  properties  of  the  thermocouple  junction 
and  the  fluid  dynamic  environment  in  which  the  element  is  used.  It  is 
important  to  note  that  the  thermocouple  output  has  to  be  properly 
compensated  for  its  time  constant  to  yield  accurate  results.  This  leads 
to  the  additional  tasks  of  determination  of  thermocouple  time  constants 
and  the  compensation  of  thermocouple  signals  for  the  time  constants. 

In  an  earlier  program'  a new  method  was  developed  to  find  the  time 
constants  under  true  operating  conditions.  Also  in  the  present  program. 
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the  compensation  is  done  later  during  data  analysis  through  suitable 
modification  of  the  Fourier  analyser  programming.  The  details  of 
this  method  are  described  in  Ref,  (14),  Neglecting  small  variations, 
the  time  constants  of  the  thermocouples  used  in  the  present  program 
are  found  to  be  5 msec. 

The  temperature  measurements  of  interest  to  this  program  are 
mean  and  fluctuating  temperatures  of  hot  gases  at  the  burner  exit 
plane.  It  is  important  to  note  that  the  entropy  noise  calculations 

' 

require  quantities  averaged  over  the  burner  exit  area.  The  conversion 
of  single  point  thermocouple  measurements  to  an  area  averaged  one 
demands  the  knowledge  of  the  spatial  distribution  of  the  temperature 

/O  \ ^ 

fluctuations.  It  has  been  found  from  an  earlier  program''  ' that  an 
accurate  determination  of  the  area  averaged  temperature  requires  a 
large  number  of  thermocouples  to  be  used  at  the  burner  exit  plane. 

Limited  by  the  availability  of  data  channels,  in  the  present  program, 
an  array  of  7 thermocouples  has  been  used.  The  schematic  arrangement 
of  the  thermocouples  is  shown  in  Fig.  5, 

Velocity  Fluctuation  Measurements  | 

I 

The  theoretical  deduction  of  various  noise  source  strengths 
requires  the  knowledge  of  impedance  coefficients  at  the  nozzle  en- 
trance plane.  The  calculation  of  these  coefficients  involve  pressure^ 


temperature  and  velocity  fluctuations.  The  pressure  and  temperature 
measurement  procedures  have  been  described  in  the  proceeding  para- 
graphs, The  velocity  fluctuations  entering  the  nozzle  were  to  be  meas- 
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ured  by  cooled  film  anemometers.  Two  such  anemometers  are  positioned 
at  the  nozzle  entrance  plane.  The  cooled  film  anemometer  is  very 
similar  to  the  hot  film  one  except  that  the  former  permits  velocity 
measurements  in  high  temperature  environments.  In  the  normal  hot 
wire  or  hot  film  anemometer  measurements,  the  temperature  of  the 
sensor  must  be  above  the  temperature  of  the  environment,  whereas 
the  cooled  film  sensor  operates  below  the  environment  temperature 
through  the  addition  of  a heat  sink.  The  basic  theory  involved  in  the 
cooled  film  anemometer  measurements  and  the  impedance  calculation 
procedures  are  described  in  Appendix  B,  As  explained  there,  these 
measurements  did  not  work  out  properly. 

Data  Acquisition 

The  pressure,  temperature  and  velocity  fluctuation  data  have 
been  obtained  through  the  instrumentation  scheme  shown  in  Fig,  6,  The 
overall  sound  pressure  levels,  registered  by  microphones,  are  read  out 
on  a BiAiel  and  Kjaer  type  2606  microphone  amplifier,  one  at  a time.  All 
the  pressure  transducers,  thermocouples  and  cooled  film  anemometer 
signals  have  been  amplified  by  an  array  of  NEFF  type  122  amplifiers 
and  recorded  on  a PR-2200,  14  channel  magnetic  tape  recorder  at  a 
tape  speed  of  30  ips.  The  r.m.s,  levels  of  the  amplified  signals  are 
indicated  by  a true  r.m.s,  voltmeter  whereas  their  visual  display  is 
available  on  the  oscilloscope  screen. 


Figure  6.  Measurement  System  Schematic. 
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V.  Data  Analysis  Procedures 

All  data  were  recorded  on  a 14  channel  IM  tape  recorder  for 
later  analysis  on  a Fourier  analysis  system.  The  primary  quantities 
needed  are  cross  spectral  data  between  all  interior  sensors  and  the 
external  microphones  and  all  interior  sensors  against  themselves. 

It  was  determined  that  the  directionality  of  the  exterior 
sound  was  higher  than  expected  on  all  runs  with  a general  downstream 
propagation  preferred.  The  microphone  at  90°  to  the  axis  was  selected 
as  the  preferred  microphone  since  its  reading  was  usually  midway  be- 
tween the  upstream  and  downstream  microphones,  Sound  power  was  cal- 
culated assuming  spherical  outward  radiation  using  the  90°  microphone 
reading. 

The  first  data  check  is  to  determine  if  there  is  interior  micro- 
phone contamination  by  hydrodynamic  noise.  This  usually  only  occurs 
at  very  low  frequency  when  there  would  ordinarily  be  very  high  coherence 
between  the  two  interior  microphones,  in  the  absence  of  contamina- 
tion. Shown  in  Fig,  7 is  a typical  coherence  obtained.  It  is  concluded 
that  there  was  negligible  contamination  of  the  interior  signals  on 
all  tests.  Consequently,  the  interior  theory  was  deemed  valid  for  use 
in  extracting  source  strengths. 

However,  this  high  coherence  between  the  two  signals  caused 

some  difficulties  in  reduction  of  the  data  at  low  frequencies.  It  is 

known  that  the  accuracy  of  multiple  coherence  estimates  degenerates 

(9) 

if  the  Inputs  are  highly  correlated.'  ' In  fact,  the  data  around 
100  Hz  were  lost  on  all  runs  because  of  this  problem,  Ihis  was  used 


IS 
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Figure  7.  Coherence  between  interior  microphones 


to  advantage  in  the  calibration  of  microphones.  The  infinite  tubes 
to  the  two  interior  microphones  carry  different  temperature  gases, 
because  of  their  locations,  and  could  be  expected  to  give  different 
readings  for  the  same  true  pressures.  The  100  Hz  data  were  used  to 
generate  a correction  factor  for  one  of  the  microphones  to  bring 
the  two  readings  to  the  same  values.  This  was  possible  because  of 
the  demonstrated  high  coherence  and  the  theoretical  demand  that  at 
this  low  frequency  there  should  be  no  spatial  structure  to  the  pres- 
sure field. 

Hie  magnitude  squared  of  the  cross  spectra  measured  were  some- 
times quite  low,  of  the  order  of  0.01  times  the  product  of  the  two 
auto  spectra  of  interest.  This  routinely  happened,  for  example,  when 
taking  cross  spectra  between  two  widely  spaced  thermocouples  or  be- 
tween a thermocouple  and  a far  field  microphone.  The  accuracy  with 
i^ich  a cross  spectrum  can  be  determined  is  determined  by  the  number 
of  samples  used  in  the  ensemble  averaging  procedure.  For  example  the 
variance  in  an  estimate  of  S^j  , o?j,i8  of  the  order  of  S^^S^^/N, 
where  N is  the  number  of  samples.  To  extract  a cross  spectral  density 
where  the  magnitude  of  the  variance  is  of  the  order  of  one  tenth  of 
the  square  of  the  spectral  density  of  interest  is 
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Consequently,  N must  be  of  the  order  of  1000  to  have  any  confidence 


1 


in  the  data  measured 


In  the  first  data  run  900  averages  were  taken  by  a combination 
of  ensemble  and  frequency  averaging,  100  ensemble  averages  were  taken 
with  a frequency  bandwidth  of  8 Hz,  Then  9 adjacent  bands  were  aver- 
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aged.  After  the  first  fully  instrumented  run,  however,  it  was  seen 
that  the  multiple  coherence  was  going  to  be  lower  than  anticipated. 
Because  the  inputs  are  highly  correlated  at  low  frequency  and  the 
accuracy  of  multiple  coherence  estimates  drops  sharply  when  this  occurs, 
all  following  runs  were  reduced  with  7000  averages,  made  by  a 1000 
ensemble  average  with  7 frequency  band  averages. 

Compensation  had  to  be  applied  to  the  thermocouples  because 
their  response  starts  rolling  off  at  6 dB/octave  above  about  40  Hz, 

They  were  compensated  by  the  technique  of  Ref , (14) where  the  correc- 
tion is  applied  during  Fourier  analysis  rather  than  by  electronic 
compensation  before  recording.  The  rolloff  forces  an  upper  frequency 
limit  of  about  2500  Hz,  above  vdiich  the  electronic  noise  limit  is 
reached , 

The  basic  measurements  were  analysed  to  produce  ®x1’ 

S ,S  ,S  ,S,,S  ,S  , and  S , However,  the  measurements 
involving  t must  be  discussed  separately  since  t is  a cross  section 
average  quantity  while  multiple  point  measurements  were  made  of  tem- 
perature, In  fact,  what  was  measured  and  analysed  were  multiple  point 
cross  spectra  between  individual  thermocouples  and  the  other  measure- 
ments, Consider  the  construction  of  from  Eq,  (16),  It  becomes 
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K ^ •*  t(x)  t(jr) 


The  transformation  from  x,  ^ to  x,  T|  coordinates  with  ^ = x + T] 
and  then  shifting  to  the  cylindrical  coordinate  systems  of  Fig,  g 


yields 
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This  complicated  formula  is  simplified  first  by  assuming  the  tempera- 
ture cross  spectrum  and  mean  temperature  is  isotropic  in  0 and  cp  so 
they  are  functions  of  p and  r alone.  This  yields 
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niese  integrals  are  then  discretized  and  numerically  evaluated  from 
the  data  by  a very  coarse  integration  scheme.  In  practice,  a pair  of 
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data  at  the  separation  distances  p = 0,  1,  2 and  3 inches  is  taken. 

This  is  considered  a universal  curve  when  normalized  by  the  auto  ! 

spectrum  at  p * 0.  Consequently,  the  r value  is  meaningless  and  the 
r variation  has  been  averaged  out.  The  above  equations  can  be  re- 
written in  the  form 


where  the  ratio  A /A  , is  a rough  measure  of  the  area  fraction 

cor  e * 

over  which  the  temperatures  are  correlated.  This  computation  proce- 
dure is  illustrated  in  Fig.  9 for  a particular  run.  Striking  in  this 
figure  Is  the  fact  that  the  space-separated  temperatures  become  more 
hi^ly  correlated  at  high  frequency  than  at  low  frequency. 

Because  of  its  coarse  construction,  is  the  quantity  known 
with  the  least  accuracy.  By  construction,  since  it  is  an  auto  spec- 
trum, it  must  be  real.  However,  at  low  frequency,  two  frequency  points 

were  encountered  where  A /A  were  negative,  which  is  impossible. 

cor  e 

The  magnitude  of  this  negative  amount  was  about  .03,  giving  an  esti- 
mate of  the  resolution  of  the  method.  When  this  occurred,  A /A 

cor  e 

was  set  equal  to  zero. 

Construction  of  the  cross  spectrum  of  t with  pressure  intro- 
duced no  special  problem.  Two  points  of  temperatures  were  used  and 
the  cross  spectra  were  averaged. 

From  these  data,  the  spectra  S^,  S^,  S^, 

are  constructed  by  manipulations  on  Eqs.  (12).  S^,  , and 
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tfmpfrature  correlation  area 


Figure  9»  Typical  correlation  area  curve 
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needed  in  Eqs.  (18),  are  also  generated  by  computations  from  the 
fundamental  data.  In  the  theoretical  computations  the  a , b and  6 

^ OD  ’ (U  W 

values  are  taken  as  generated  In  Section  VI.  Ordinary  coherence 
functions  between  all  quantities  and  the  multiple  coherence  function 
are  then  computed.  The  data  are  checked  so  that  all  coherence  func- 
tions at  each  frequency  must  lie  between  0 and  1 or  else  that  fre- 
quency point  is  rejected.  Data  rejection  at  100  Hz  always  occurred.  In 
the  first  fully  instrumented  run,  some  data  rejection  at  200  and  300  Hz 
also  occurred. 

A complete  statistical  analysis  of  errors  in  the  measurements 
and  data  reduction  process  was  beyond  the  scope  of  the  program.  How- 
ever, judging  from  data  scatter  in  the  spectral  results  to  follow  and 
knowledge  of  the  accuracy  of  determination  of  some  of  the  cross  spec- 
tral properties  leads  one  to  the  following  judgements:  a)  sound  power 
results  are  probably  correct  to  +1.5  dB,  b)  spectral  determination 
of  source  quantities  is  probably  accurate  to  + 3 dB  and  c)  mean  square 
determination  of  source  quantities  is  probably  correct  to  + 2 dB. 
Further  discussions  of  the  split  of  the  sound  power  into  various  caus- 
al elements  is  deferred  until  later. 
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VI.  Admittance  and  Transmittance  Determinations 
It  has  been  mentioned  in  Section  II  that  application  of  theory 


to  the  data  reduction  process  requires  knowledge  of  the  admittance 

(3) 

coefficients  a , b and  ft  » In  prior  calculations  and  tests  the 

UJ*  (0 

effect  of  ft  was  found  to  be  small,  as  long  as  it  is  small.  It  was 

found  8^  0,06  was  adequate  for  matching  of  theory  and  data  and 

this  value  will  be  used  here.  However,  a and  b have  a profound  in- 

U)  U) 

fluence  on  the  deduction  of  the  vorticity  noise  strength  and  they 

must  therefore  be  known  with  some  precision,  Tliree  methods  were  used; 

these  were  a)  theoretical  calculations  using  Candel's  theory 

b)  measurements  using  a combination  of  cooled  film,  thermocouple  and 

pressure  measurements  and  c)experiments  on  a model  nozzle  in  an  entropy 

tube  facility.  These  will  be  described  below. 

The  diffuser  configuration  employed  was  to  have  accomplished 

two  things;  these  are  a)  to  drop  the  jet  noise  contamination,  which 

it  did  and  b)  to  act  as  an  impedance  matching  horn,  which  it  only 

partially  did.  The  hope  was  that  all  acoustic  energy  passing  the 

choked  throat  would  get  out  to  the  far  field.  Apparently,  it  did  not 

and  the  transmission  properties  of  the  nozzle  diffuser  assembly  are  of 

interest  since  they  influence  the  results.  Moreover,  the  transmittances 

due  to  entropy  and  vorticity  incident  on  the  nozzle  cannot  be  explained 

through  a and  b alone  because  there  is  continuous  generation  of  re- 
U)  U) 

fleeted  and  transmitted  waves  as  the  incident  entropy  and  vorticity 
traverse  the  nozzle.  Consequently,  the  transmittance  due  to  incident 
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Isentroplc  and  entropy  waves  was  looked  at  in  some  detail  as  an  aid 
in  interpretation  of  results. 

The  theoretical  calculations  of  admittance  coefficients  fol- 
lows from  Ref.  (15)  where  the  mean  velocity  is  assumed  linear  in 
distance  from  the  entrance  plane  to  the  throat.  For  one  of  the  nozzles 
a theoretical  calculation  is  shown  in  Fig. 10,  The  values  undergo 
a monotonic  decline  in  absolute  magnitude  with  an  increase  in  frequency. 
The  a^  magnitude  undergoes  a rise  with  frequency  in  the  frequency  range 
shown.  Actually,  as  the  wavelength  becomes  comparable  to  the  nozzle 
length,  |a^j  would  undergo  a maximum.  However,  the  frequencies  con- 
sidered in  this  program  are  low  enough  that  the  behavior  of  Fig.  10 
is  all  that  need  be  considered.  The  basic  point  to  be  noted  is  that 

la  ! and  lb  , are  small  nxunbers,  of  the  order  of  the  mean  entrance 
• u)‘  * U)' 

Mach  number.  Ihese  low  admittance  values  mean  that  the  nozzle  is 
highly  reflective  to  incident  isentropic  waves  but  quite  sensitive 
to  incident  entropy  waves,  insofar  as  creation  of  pressure  waves  is 
concerned.  On  the  other  hand,  as  frequency  increases  the  nozzle  be- 
comes a)  less  reflective  to  incident  isentropic  waves  and  b)  less 
sensitive  to  generation  of  pressure  by  incident  entropy  waves. 

Confirmation  of  the  theoretically  expected  behavior  was  given 
by  the  results  of  the  entropy  tube  tests  which  are  fully  described 
in  Appendix  A.  A model  nozzle,  designed  for  the  same  axial  Mach  num- 
ber gradient  as  the  hot  flow  # 2 nozzle,  was  tested  in  the  entropy 
tube  and  bombarded  with  isentropic  and  entropy  waves.  Unfortunately 
the  signal  to  noise  ratio  was  not  good  at  frequencies  above  250-500 
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Figure  lO.  Isentropic  and  non-isentropic  admittance 
coefficient  calculated  theoretically  for 
, particular  nozzle. 


Hz  so  that  the  results  are  relegated  to  Appendix  A.  Reference  to 
Figs,  14  and  22  of  that  Appendix  show  the  expected  trends  but 
quantitative  comparison  is  not  highly  satisfactory. 

The  actual  hot  flow  test  nozzle  passes  a flow  which  is  not 
uniform  in  temperature  due  to  a bypass  cooling  flow  which  is  at  the 
wall#  Consequently,  it  was  desired  to  measure  the  admittance  coeffi- 
cients directly  in  the  hot  flow  case.  In  Appendix  B are  described 
the  results  of  that  effort.  Unfortunately  the  claims  of  the  manu- 
facturer concerning  a piece  of  instrumentation  (cooled  film  probe) 
were  shown  to  be  highly  inaccurate.  As  a consequence,  the  attempted 
measurements  failed. 

Because  of  the  above  problems  it  was  decided  to  use  the  theo- 
retical admittance  values  in  the  data  reduction  procedure.  The  justi- 
fication is  through  the  trends  exhibited  in  the  entropy  tube  results, 
the  fact  that  when  la  | and  |b  | are  low  enough  they  affect  the  re- 

' U)'  ' U)i 

suits  in  only  a minor  way,  and  the  fact  that  nozzles  have  behaved 
well  when  compared  against  theory  in  other  programs . 

Concerning  the  transmittance  properties  of  the  nozzle-diffuser 
combination.  Figures  18  and  23  in  Appendix  A give  a feel  for  what  may 
be  expected.  As  frequency  rises  and  the  isentropic  admittance  increases, 
isentropic  waves  having  an  easier  time  getting  through.  Thus,  it  is 
expected  that  the  transmittance  will  rise  with  frequency.  More  com- 
plex behavior  with  respect  to  the  entropy  transmittance  is  expected, 
however,  as  the  entropy  admittance  coefficient  drops  with  a frequency 
rise,  the  pressure  generation  ability  of  the  entropy  wave  drops.  On 
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the  other  hand,  the  transmittance  of  the  generated  pressure  wave  in- 
creases. Figure  18  of  Appendix  A confirms  that  there  are  competing 
mechanisms. 

The  actual  transmittance  coefficients  in  the  hot  firing  were 
obtained  by  comparison  of  external  and  internal  pressures.  First, 
the  theory  of  Eqs.  (12)  has  to  be  reworked  to  define  the  incident 


pressure  pulse,  since  transmittances  are  defined  in  terms  of  incident 
and  transmitted  pulses.  Hie  theory  as  presented  in  Eqs.  (12)  gives 
the  complete  pressure  at  a point,  consisting  of  incident  and  reflected 
pulses.  Decomposing  the  solution  of  Eqs.  (12)  into  left  and  right 
running  waves,  the  wave  incident  on  the  nozzle  is  given  by 


ji 


results.  Here  H. , is  measured  and  H,  is  theoretically  computed. 

A1  Ao 

For  the  entropy  transmittance  coefficient  the  dimensionless  definition 
PiUi  /p(®) 

- A __| — ...  ^ 


I-Wl  • K •*  ^ 

is  used  ^ , where  one  has  to  be  careful  here  to  ex- 

U) 

elude  from  pu^j^  components  due  to  the  isentropic  incident  wave  caused 
by  reflections  of  entropy  generated  isentropic  waves.  Doing  this, 


T.-  + H , 

il  ctI 


The  quantities  of  Eqs.  (24)  and  (25)  are  also  generated  during 
the  data  reduction  process  and  are  shown  in  Figs.  11  and  12, 

shows  the  striking  rise  in  transmittance  with  a frequency  rise, 
as  expected,  but  it  does  this  for  both  the  jet  and  diffuser  termina- 
tions. The  second  feature  to  note  is  that  the  jet  transmittance  is 
substantially  lower  than  for  the  three  diffusers.  Consequently,  there 
is  some  mechanism,  which  is  not  understood  at  this  time,  which  favors 
sound  transmittance  in  diffused  flows  following  a choking  point.  This 
will  show  up  later  as  a sound  level  significantly  higher  than  can  be 
explained  by  the  correlation  efforts  of  Section  III,  The  increase  in 
transmittance  with  frequency  is  characteristic  of  an  impedance  match- 
ing horn,  but  the  jet  effect  is  not  understood. 

As  expected  from  Appendix  A,  Fig,  12  shows  no  particularly  strik- 
ing behavior  for  the  entropy  transmission  coefficient.  The  jet  again 
shows  a depressed  transmittance  curve  with  a shape  which  is  similar 
to  the  nozzle-diffuser  curves. 

In  Figs,  11  and  12  only  magnitude  has  been  presented  since 
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ISENTRa  . TRANSMISSION  COEFFOENTS 


FREQUENCY  (HZ) 


ENTROPY  TRANSMISSION  COEFFIQENTS 


Entropy  transmission  coefficients  for  the 
four  fully  reduced  runs. 
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phase  information  depends  upon  the  distance  between  the  source  and 
receiver  points.  The  vertical  scale  also  should  be  viewed  on  a rela- 
tive magnitude  basis  since  the  absolute  magnitude  depends  upon  the 
location  of  the  exterior  microphone.  Only  the  shape  on  the  frequency 
plot  is  of  interest  here. 


I 


VII.  Experimental  Results 
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Source  Identification 

The  spectral  results  to  be  presented  are  two-sided  spectra  cor- 
rected to  a 100  Hz  bandwidth.  Consequently,  the  mean  square  value  of 
a quantity  may  be  obtained  by  band  s\immation,  and  then  multiplication 
by  2,  Moreover,  the  spectra  are  always  of  normalized  quantities  so 
that  mean  square  values  are  really  fractional  mean  square  quantities. 

Shown  in  Fig.  13  is  a typical  auto  spectrum  of  an  interior  mi- 
crophone and  in  Fig,  14  is  shown  a typical  compensated  thermocouple 
spectrum.  The  microphone  resonances  correspond  to  the  nearly-closed- 
end  longitudinal  resonances  in  the  vicinity  of  800,  1600,  etc.  Hz, 

In  addition  there  is  a peak  near  100-200  Hz  which  has  a counterpart 
in  the  thermocouple  spectrum;  this  is  perhaps  the  entropy  wave- 
acoustic  wave  resonance  which  was  shown  to  exist  in  Ref,  (3),  Tlie 
thermocouple  spectrum  is  not  "closed"  in  the  sense  that  there  is  high 
confidence  that  the  area  under  the  curve  will  give  the  mean  square 
value  of  the  relative  temperature  fluctuation.  However,  over  the  range 
of  frequencies  shown  the  rms  fluctuation  is  of  the  order  of  12-13%  of 
the  mean  temperature.  Recall,  the  reason  higher  frequencies  are  not 
shown  is  that  above  2500  Hz  the  thermocouple  output  is  in  the  elec- 
tronic noise.  This  will  not,  however,  affect  any  of  the  conclusions 
of  the  work.  In  Fig,  13,  the  rapid  falloff  with  frequency  is  expected 
to  yield  a rapidly  diminishing  combustion  noise  source  strength,  in 
accordance  with  Eq,  (15),  This  will  indeed  be  found  to  be  the  case. 


PRESSURE  FLUCTUATION  SPECTRUM 


Figure  13,  Auto  spectrum  of  an  interior  microphone 
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pectrum  of  an  exit  plane  thermocoup li 


In  Fig,  15  are  displayed  the  auto  spectra  of  the  90°  exterior 
microphone  for  the  four  basic  runs.  The  extent  and  detail  of  the  high 
frequency  information  increased  with  the  order  in  which  the  runs  were 
made  for  reasons  explained  later.  For  the  three  runs  with  diffusers 
there  is  a general  level  increase  and  shift  to  higher  frequencies  as 
the  mass  flow  is  increased.  The  most  striking  effect  is  the  depression 
of  the  noise  when  there  was  no  diffuser.  The  jet-only  spectrum  is  still 
climbing  at  2500  Hz  because  of  jet  noise  masking  of  the  combustion 
noise  results,  as  will  be  seen.  The  spectra  of  Fig.  15  contains  noise 
from  all  possible  sources,  not  just  those  inside  of  the  combustor; 

Fig.  15  is  lerely  raw  noise  data  and  must  be  corrected  to  obtain  noise 
caused  by  Interior  sources. 

Shown  in  Fig,  16  are  the  correlation  areas  for  the  temperature 
fluctuations.  Recall  these  are  needed  to  compute  the  area-average  of 
the  temperature  fluctuations.  There  were  two  problem  points  - at  300 
Hz  on  the  first  run  and  at  400  Hz  on  the  last  run  the  correlation  area 
was  negative,  which  is  a physical  and  mathematical  impossibility  and 
reflects  the  use  of  a coarse  net  of  thermocouples.  These  problem  points 
will  cause  spurious  results  to  follow,  but  overall  will  cause  little 
effect  since  they  are  in  but  one  band  of  the  spectrum.  There  are  two 
interesting  results  from  Fig,  16,  These  are  a)  the  correlation  area 
"resonates"  at  the  can  acoustic  mode  frequencies  and  b)  the  correla- 
tion area  generally  increases  with  frequency.  The  resonance  was  ex- 
pected since  the  temperature  fluctuation  is  made  up  of  two  parts  - 
one  due  to  random  hot  spots  and  the  second  due  to  the  wave  motion  in 
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the  can.  The  wave  motion,  being  basically  a plane  wave,  is  correlated 
over  the  entire  can  diameter  and  dominates  the  temperature  fluctua- 
tion at  resonance.  The  general  rise  with  frequency,  however,  is  sur- 
prising. One  generally  associates  higher  frequencies  with  smaller 
hot  spots.  This  result  suggests  that  although  the  longitudinal  length 
scale  may  decrease  with  an  increase  in  frequency  the  lateral  scale 
increases.  Within  the  limits  of  data  scatter  expected,  all  curves  of 
Fig,  16  are  believed  consistent.  There  was  no  difference  expected 
between  the  runs  for  the  correlation  area  with  only  a factor  of  2:1 
variation  in  mass  flow. 

Shown  in  Fig,  17  are  the  multiple  coherence  functions  for  the 
four  basic  runs.  It  is  this  function  which  shows  the  fraction  of 
noise  on  Fig,  15  which  is  linearly,  causally  related  to  the  interior 
sources.  There  is  a general  rise  in  coherence  with  frequency,  except 
for  the  nozzle-only  case  which  had  heavy  jet  noise  masking  at  high 
frequency.  In  any  event  at  low  frequency  there  is  always  relatively 
poor  multiple  coherence  indicating  a)  nonlinear  behavior,  b)  measure- 
ment noise  or  c)  the  presence  of  other,  unaccounted  for,  sources. 
However,  the  multiple  coherence  level  is  of  such  magnitude  that  there 
is  statistical  confidence  in  the  results,  even  at  low  frequency.  While 
it  only  has  meaning  as  an  indicator  of  average  coherence,  the  band 
summation  averages  of  the  multiple  coherence  are  0,28,  0.35,  0.14  and 
0.27,  in  the  order  of  increasing  mass  flow. 

The  coherent  noise  spectrum  is  obtained  by  multiplying  the  re- 
sults of  Fig,  17  by  the  results  of  Fig,  15.  The  results  are  shown  in 
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MULTPLE  COHERENCE 


Figure  17.  Multiple  coherence  function  for  the  four 
basic  runs. 


Figure  18.  Coherent  noise  spectra  of  the  90  microphone 
for  the  four  basic  runs. 


Fig,  18.  The  reason  for  the  ever-widening  frequency  range  of  the  analy- 
sis is  now  apparent.  In  the  first  run  (0,18  kg/sec)  the  coherent  noise 
spectrum  is  not  closed  in  the  sense  that  the  spectrum  is  still  rising 
at  the  upper  analysis  limit  of  1500  Hz.  From  all  prior  programs  it  was 
anticipated  that  the  coherent  noise  would  be  well  down  from  the  maxi- 
mum by  about  1000  Hz,  This  did  not  occur  and  forms  the  basis  for  a 
major  result  of  the  program.  The  analysis  was  widened  to  2000  Hz  for 
the  second  run  (no  diffuser)  but  still  the  spectrum  did  not  close. 

For  the  last  two  runs  the  analysis  was  widened  to  4000  Hz  but  it  was 
found  that  above  2500  Hz  there  was  a rapid  drop  in  the  coherent  spec- 
trum so  that  results  are  only  presented  to  2500  Hz,  Most  of  the  co- 
herent noise  is  contained  in  the  "hianp"  around  2000  Hz  and  this  is  a 
new  finding  for  combustor  noise. 

The  data  below  500  Hz  for  the  coherent  noise  are  somewhat  sus- 
pect because  of  the  low  multiple  coherence  and  problems  with  the  tem- 
perature correlation  area.  However,  the  low  frequency  range  is  not 
the  major  noise  contributor  so  that  this  problem  is  of  little  conse- 
quence, There  is  a general  overall  increase  in  level  as  mass  flow 
increases,  as  expected  and  the  results  of  Fig.  15  persist  that  the 
jet-only  run  is  at  a much  lower  level  than  the  runs  with  the  diffuser. 
If  the  basic  source  strength  within  the  combustor  remain  relatively 
fixed,  as  expected,  this  level  change  is  consistent  with  the  trans- 
mittance results  of  Figs.  11  and  12. 

Consider  now  the  derived  source  strengths.  These  are  the  auto 

spectra  S.,,  S and  S . They  are  presented  in  Figs.  19-21.  S,, 
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ENTROPY  SOURCE  STRENGTH 


20.  Entropy  noise  source  spectrum  for  the  four 
basic  runs. 


VORTIOTY  SOURCE  STRENGTH 
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Figure  21.  Vorticity  noise  source  spectrum  for  the 
four  basic  runs. 
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has  a maximum  at  low  frequency,  rapidly  falls  then  has  a hump  around 

2000  Hz,  The  major  contribution  to  the  mean  square  value  of  A comes 

(3) 

at  low  frequency.  In  a prior  program  the  derived  peaked  at 
200  Hz,  but  recall  here  that  the  data  at  100  Hz  were  lost  so  that 
this  fact  could  not  be  confirmed  here.  The  same  general  trend  holds 
for  all  curves,  as  expected,  with  the  only  major  discrepancy  being 
a hump  in  the  high  mass  flow  spectrum  around  1000  Hz, 

In  Fig,  20  the  entropy  source  strength,  S , is  displayed  and 

CTO 

there  is  high  consistency  between  runs,  as  expected.  In  Fig,  21, 
shows  the  same  general  trends  between  the  four  runs,  but  there  is  an 
increase  in  overall  level  with  an  increase  in  mass  flow.  While  S 

CTCT 

is  relatively  flat,  shows  roughly  the  same  behavior  as  the 

curves  of  Fig,  19, 

The  overall  behavior  of  the  mean  square  values  of  the  source 
strength  is  obtained  by  band  summation  and  is  shown  in  Table  3, 


Table  3 

Mean  Square  Values  of  the  Interior 
Source  Strengths 


Run 

A 

CM 

«< 

V 

2 

<CT  > 

A 

CM 

> 

V 

0,18  kg/sec  - diffuser 

1,73  X lO"^ 

8,26  X lO"^ 

1.45  X lo"^ 

0,25  kg/sec  - diffuser 

1,02  X lO"^ 

1,65  X lo'^ 

1.01  X 10“^ 

0,34  kg/sec  - diffuser 

1.80  X lo"^ 

2,98  X 10“^ 

1.76  X lo"^ 

In  Table  3 there  is  a clear  trend  for  the  mean  square  values  of  the 
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source  strengths  to  rise  with  mass  flow,  which  was  not:  expected  in 

2 2 

the  theory  of  Section  II,  As  an  example  < a > “ m , This  kind  of  a 

2 2 

trend  is  apparent  in  the  bottom  two  entries  for  < A > and  < V >, 

However,  the  first  entry  for  these  two  quantities  indicates  a stronger 

2 2 

dependence.  On  the  other  hand,  < A > and  < V > are  highly  influenced 

by  the  low  frequency  portions  of  the  spectrum  (see  Figs,  19  and  21) 

where  the  data  is  least  accurate  and  where  there  is  little  coherent 

noise  output.  At  higher  frequencies  there  is  less  of  a difference 

between  the  first  and  last  runs.  Because  of  the  relative  insensitiv- 

2 

ity  of  S to  frequency  the  deduction  of  the  A dependence  is  pre- 
aa 

ferred  for  all  source  strengths.  This  would  bring  the  theory  in  closer 
line  with  the  empirical  correlation  of  Section  III,  although  the  physics 
of  why  this  trend  occurs  is  unclear.  At  any  rate,  this  interpretation 
will  be  accepted  and  the  conjecture  of  Section  III  that  the  additional 
dependence  of  power  on  as  compared  with  theory,  is  due  to  aero- 

dynamic noise  contamination  will  be  dropped. 

Now  the  question  arises  as  to  the  separability  of  the  coherent 
noise  into  components  due  to  the  various  sources.  The  first  bit  of 
damning  evidence  comes  from  Fig,  22,  where  the  ordinary  coherence  be- 
tween sources  is  calculated.  First  of  all,  at  low  frequencies  all 
sources  are  highly  coherent  with  each  other.  This  was  expected  between 
entropy  and  combustion  noise  from  Ref,  (3),  but  it  is  now  shown  to  in- 
clude vorticity  as  well.  At  the  higher  frequencies  (>2000  Hz)  there  is 
an  exceedingly  high  coherence  between  the  combustion  and  vorticity 
sources.  In  the  mid  range  there  is  often  a high  coherence  between 


COHERENCE  ESTIMATES 


Figure  22.  Ordinary  coherence  between  sources  for  the 
four  basic  runs. 
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vorticity  and  entropy.  These  facts  are  merely  stating  that  the  off 
diagonal  elements  of  Eq.  (18)  are  important  to  the  overall  sound 
that  is  created,  and  one  cannot  justifiably  split  the  sound  into 
unique  components. 

As  an  example  consider  that  only  two  sources  are  present,  both 
of  which  are  in  phase  with  the  radiated  sound  if  each  were  to  act 
alone.  If  they  are  correlated  and  out  of  phase  with  each  other  there 
would  be  no  radiated  sound,  but  if  they  were  in  phase  there  would  be 
twice  as  much  sound  as  compared  with  the  sum  of  the  sound  from  each 
individual  source.  Given,  however,  that  there  is  a strong  desire  to 
say  \diich  cause  is  dominant  in  producing  sound  and  that  the  coherence 
between  sources  is  rarely  perfect  on  Fig,  22,  an  approximate  procedure 
will  be  used  to  split  the  sound. 

The  ordinary  coherence  between  the  sources  and  the  measured 
sound  was  computed  and  multiplied  by  the  measured  sound,  and  these 
values  are  called  the  "sound  due  to  source  i",  whatever  source  i is. 

The  split  is  then  displayed  on  Figs,  23-26,  Ignore  for  the  moment  the 
bold  square  points  on  Figs,  23  and  24,  If  one  adds  up  the  components 
one  will  sometimes  get  greater  than  and  sometimes  less  than  the  upper 
curve  which  is  the  multiple  coherent  noise.  This  is  due  to  the  prob- 
lem above  that  the  sources  are  coherent  with  each  other.  Nevertheless, 
approximately,  two  major  facts  are  seen.  These  are  a)  at  high  frequency 
(say  >1100  Hz)  vorticity  and  combustion  noise  are  dominant  and,  recall, 
usually  highly  correlated  and  b)  the  only  frequency  regime  where  entropy 
noise  is  strong  is  at  low  frequency  (say  < 600  Hz). 
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Figure  25,  Multiple  coherent  noise  and  the  split  of  the 
noise  into  the  various  causal  elements  for 
the  0.30  kg/sec  run. 


The  comparisons  become  much  less  ragged  if  one  does  a band 
summation  on  the  various  components.  If  one  adds  the  bands  of  the 
various  components  and  divides  by  the  summation  of  the  bands  of  the 
overall  multiple  coherence  runs,  the  figures  0.90,  0.65  and  0.91  re- 
sult for  the  three  diffuser  runs  in  order  of  increasing  mass  flow. 

That  is,  the  components  add  up  to  less  than  the  whole,  but  comparable 
to  it,  because  of  the  coherence  between  sources. 

Shown  in  Table  4 are  summations  of  the  various  components  from 
200-1500  Hz  and  200-2500  Hz  expressed  as  a fraction  of  the  total  of 
all  components.  The  1500  Hz  point  is  selected  because  the  first  run 
only  was  analysed  up  to  1500  Hz.  In  the  last  two  runs  it  is  clear 
that  combustion  and  vorticity  noise  dominate  above  1500  Hz  and,  in 
fact,  the  second  run  had  this  domination  at  1500  Hz.  The  first  run 
had  roughly  equivalent  source  causality  at  1500  Hz.  Based  on  the  last  j 

j 

two  runs,  however,  and  the  apparent  takeover  of  combustion  and  vor- 
ticity at  1500  Hz  on  the  first  run  (see  Fig.  23),  it  can  be  surmised 
that  if  the  first  run  had  been  analysed  out  . 2500  Hz,  the  same  split 

would  have  occurred  as  with  the  last  two  runs.  Such  figures  are  en- 
tered in  parentheses  in  Table  4 as  an  average  from  the  last  two  runs.  . ! 

The  major  conclusions  here  are  that  a)  combustion  and  vorticity  noise 
dominate  the  noise  output  and  are  highly  correlated  and  b)  only  at  low 

frequency  is  entropy  noise  a strong  source.  What  is  being  said  here  is  \ 

that  iji  this  apparatus  the  presence  of  vorticity-uozzle  interaction  adds 
about  a factor  of  two  to  the  overall  noise  as  compared  to  vAiat  would 
occur  if  only  combustion  noise  were  present.  The  relative  split  be- 


\ 

c 


tween  the  sources  in  the  nozzle-only  run  was  63%  , 23%.  aad  14%  and 
the  first  two  would  undoubtedly  have  climbed  relative  to  the  last 


one  if  the  analysis  had 

been  carried 

to  2500  Hz. 

Table  4 

Split 

of  the  Noise 

into  Causal 

Properties  for  the  Diffuser  Runs 

(Hz) 

% 

Run 

Freq. 

Combustion 

Vorticity 

Entropy 

0,18  kg/sec 

1500 

24 

40 

36 

2500 

(45) 

(51) 

(4) 

0.25  kg/sec 

1500 

40 

53 

7 

2500 

42 

56 

2 

0*34  kg/f«c 

1500 

32 

43 

25 

2500 

49 

46 

5 

Returning 

for  a moment  to  Figs, 

23  and  24,  the 

derived  sound 

due  to  vorticity  is  displayed  by  the  bold  squares  \dien  the  impedances 

measured  by  the  cooled  film  probes  are  used.  These  are  not  believed, 

however,  for  the  reasons  cited  in  Appendix  B,  Moreover,  the  indicated 
2 

<v  > rises  by  almost  a factor  of  100  which  yields  an  unrealistic  rms 
fluctuation  level  (nearly  100%)  for  an  area-averaged  quantity. 

Sound  Power  Scaling 

The  multiple  coherent  sound  power  and  the  sound  power  scaling 
laws  are  obtained  by  multiple  corrections  to  the  data.  The  power  de- 
sired is  the  effective  sound  power  so  that  it  may  be  compared  with  the 
prediction  formula  of  Section  III, 


Consider  first  the  sound  power  for  the  four  primary  runs  where 

the  multiple  coherent  power  is  to  be  placed  on  an  effective  basis. 

Band  summation  was  applied  to  the  coherent  mean  square  pressure  of 

the  90°  microphone,  but  the  spectra  were  corrected  by  a correction 

curve  closely  related  to  Fig.  1.  It  was  assumed  that  combustion 

2 2 

noise  was  dominant  and  the  was  compared  with  for  a 

reflection  free  termination.  The  coherent  p^  was  assumed  proportional 
to  p^ , as  it  would  be  if  only  combustion  noise  was  present.  The  re- 
sult is  that  7.9  dB  must  be  added  to  the  measured  mean  square  coherent 
< p^  > to  bring  it  to  an  equivalent  reflection  free  situation.  This 

magnitude  may  have  been  expected  from  Fig.  1.  It  was  then  assumed 

2 

that  a monopole  radiation  pattern  existed  with  < p^  > then  converted 
to  sound  power.  However,  a further  correction  had  to  be  applied  to 
the  data  of  the  first  run  since  it  was  analysed  only  out  to  1500  Hz. 

It  was  noted  in  the  second  two  diffuser  runs  that  the  ratio  of  the 
2500  Hz  value  of  < p^  > to  the  1500  Hz  value  was  a consistent  2.3. 
Consequently,  this  factor  was  multiplied  to  the  sound  power  of  the 
first  run.  The  jet-only  run  was  left  unmodified,  since  at  2000  Hz, 
its  upper  analysis  limit,  the  other  runs  had  started  to  decay  in 
the  coherent  In  all  cases,  the  resultant  numbers  were  multiplied 

by  0.45,  0.42,  0.49  and  0.63,  respectively,  to  convert  the  total 
sound  to  combustion  noise  sound,  because  the  prediction  formula  is 
only  good  for  pure  combustion  noise. 

The  results  for  the  basic  runs  are  shown  in  Table  5 and  are  quite 
interesting.  The  measured  effective  coherent  sound  p)ower  is  about  10  dB 
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higher  than  the  prediction  formula  for  the  runs  with  the  diffuser, 
although  the  scaling  with  mass  flow  appears  as  predicted o On  the 
other  hand,  the  nozzle-only  run  with  a high  velocity  exit  jet  shows 
a good  comparison  with  the  formula  value.  It  is  clear  that  the  dif- 
fuser is  doing  something,  but  exactly  what  is  not  known.  The  first 
run  value  appears  slightly  high  in  comparison  with  the  other  two 
diffuser  runs,  but  it  should  be  remembered  that  the  diffuser  flow 
was  separated.  Consequently,  there  may  have  been  some  difference  in 
sound  transmission  behavior  or  additional  coherent  sound  generation. 

Table  5 

Sound  Power  from  Experiment  and  from  the 
Correlation  Formula  for  the  Four  Basic  Runs 


Acoustic  Power  (dB 

re  10  ^^w) 

Run 

Correlation  Formula 

Experiment 

0,18  kg/sec  diffuser 

111.6 

122.6 

0.25  kg/sec  diffuser 

116.4 

125.0 

0.34  kg/sec  diffuser 

120.9 

128.2 

0,30  kg/sec  nozzle  only 

116.1 

113.4 

For  the  mid-contraction  ratio  case  a set  of  runs  were  made  with 
the  diffuser  while  varying  the  airflow  and  fuel/air  ratio  off  of  the 
nominal  operating  condition.  The  90°  B&K  microphone  was  used  as  the 
reference  microphone,  monopole  radiation  was  assumed  and  the  "average" 
baseline  multiple  coherence  value  of  0.35  was  applied  directly  to  the 
overall  < p^  >.  Moreover,  the  value  from  Table  4 of  42%  combustion 
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noise  was  applied  to  the  value  and  the  7.9  dB  correction  factor  was 
also  applied.  The  results  are  shown  in  Table  6,  Notice  the  close 
agreement  between  the  band  summation  value  in  Table  5 and  the  base- 
line pKJwer  obtained  by  manipulation  of  the  overall  measured  pressure 
in  Table  6.  Uie  fuel/air  ratio  scaling  was  as  expected.  However,  the 
observed  dependence  on  mass  flow  is  stronger  than  predicted.  This 
may  be  giving  a clue  that  the  diffuser  is  indeed  at  fault  in  causing 
higher  overall  levels  than  expected.  As  mass  flow  increases  from  the 
baseline  value  a shock  must  form  in  the  diffuser  which  may  be  sensi- 
tive to  upstream  disturbances  in  causing  more  noise.  On  the  other 
hand,  as  the  mass  flow  drops  from  the  choking  condition  the  experi- 
mental sound  power  is  dropping  rapidly  toward  the  formula  value. 
However,  if  this  were  to  be  a valid  explanation  then  the  same  ex- 
treme sensitivity  would  have  had  to  have  been  present  when  fuel/air 
ratio  was  varied,  and  it  was  not. 

Table  6 

Sound  Power  from  Experiment  and  the 
Correlation  Formula  for  the  Variations 
about  the  Baseline  Case  for  Nozzle-Diffuser 
Assembly  #2 

Power  (dB  re  10  w) 


Airflow  (kg/sec) 

Fuel /Air 

Correlation 

Measured 

0.25 

0.0180 

116.4 

125.3 

0.21 

0.0175 

114.2 

119.7 

0.31 

0.0173 

118.3 

129.8 

0.25 

0.0156 

115.9 

125.2 

0.25 

0.0203 

117.1 

126.2 
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Discussion 

Upon  review  of  the  data  the  following  summary  is  offered: 

a)  The  source  quantities  behaved  similarly  in  all  runs,  except 
for  an  empirical  discovery  that  the  level  is  mass  flow  sen- 
sitive which  only  serves  to  strengthen  the  agreement  be- 
tween the  interior  theory  and  experiment. 

b)  The  impedance  of  a choked  nozzle  tends  to  favor  higher  fre- 
quencies in  the  radiated  sound  than  heretofore  observed. 

c)  There  is  high  correlation  between  the  vorticity  and  combus- 
tion noise  sources  and  they  are  dominant  in  producing  co- 
herent noise,  especially  at  the  higher  frequencies. 

d)  Hie  multiple  coherence  for  the  nozzle-diffuser  assemblies 
was  high  at  high  frequencies  where  the  majority  of  sound 
radiation  occurred. 

e)  The  coherent  noise  from  the  nozzle-only  configuration  was 
somewhat  lower  than  predicted  but  the  coherent  noise  from 
the  nozzle-diffuser  assemblies  was  much  higher  than  pre- 
dicted. 


f)  The  split  between  vorticity  and  combustion  noise  favored 
vorticity  noise  more  in  the  case  of  the  nozzle-diffuser 
apparatus  as  compared  with  the  nozzle-only  case. 


AO-A067  993 

UNCLASSIFIED 


GEORGIA  INST  OF  TECH  ATLANTA  SCHOOL  OF  AEROSPACE  ENG — ETC  F/G  20/1 

combustion  noise  prediction:  entropy/vorticity  effects. (U) 

MAR  79  W c STRAHLE*  M MUTHUKRIShNAN  OOT-FA77WA-4077 

FAA-RD-79-1  ML 


2 OF  2 

AO 

A067993 

■ 

- - 

i 

- 

1 

i ; 

- 

.. 

AL-  * 

— 

,! 

1 

- ■ ■ , 

■■  ; 

- 

1 > 

1 ■' 

1— - 



:x 

'yivi' 

END 

DATE 

FILMED 

6-79 

DOC 

was  carried  out  to  only  2000  Hz  and  there  would  have  been  some  sound 
in  the  2000-2500  Hz  band  to  have  added  to  the  overall  level.  Secondly, 
near  field  microphones  were  used  and  the  high  velocity,  high  tempera- 
ture jet  has  a high  impedance  mismatch  with  the  surroundings.  The  jet 
may  have  been  acting  as  a waveguide  for  a few  diameters  and  the  near 
field  microphones  saw  a smaller  fraction  of  the  overall  noise  than 
would  have  been  seen  in  the  far  field.  Consequently,  there  is  believed 
no  real  discrepancy  between  expectation  and  experiment  for  the  nozzle- 
only  case. 

The  nozzle-diffuser  assembly  cases  are  more  puzzling,  and  the 
conclusion  appears  inescapable  that  addltonal  noise  was  being  created 
in  the  diffuser  by  the  vortlcity  which  was  coherent  with  the  can- 
interior  vorticity.  This  is  indicated  by  the  fact  that  vortlcity  noise 
was  more  dominant  in  the  diffuser  cases  as  compared  with  the  nozzle 
case.  Thus,  a flow  with  an  adverse  pressure  gradient  may  be  susceptable 
to  vorticity-f  low  interaction  which  causes  noise-.  A second  reason  for 
increased  noise  is  perhaps  the  high  coherence  between  combustion  and 
vortlcity  noise.  As  described  above,  two  highly  coherent  sources  can 
produce  more  noise  than  the  noise  produced  by  the  sum  of  two  uncorre- 
lated sources. 

The  lack  of  an  entropy  noise  dominance  was  surprising.  However, 
it  was  due  to  the  high  frequency  dominance  of  the  noise.  Entropy  noise 
is  as  important  as  the  other  two  types  at  lower  frequencies. 
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Vlll.  Corrections  to  the  Prediction  Formula 
Because  of  the  uncertainties  in  some  of  the  results,  the  fact 
that  the  dominance  of  relative  contributions  Is  highly  frequency  sen- 
sitive, and  the  combustor  termination  appears  to  strongly  influence 
the  noise  output,  it  appears  difficult  to  give  a general  correction 
to  the  combustor  rig  formula  for  sound  power.  However,  a form  \rfiich 
is  suggested  is 

f Vr  1 ^ \ (26) 


The  first  computation  is  that  of  from  Eq.  (21)  using  the 
empirical  coefficients  of  Table  1,  This  gives  the  power  which  would 
be  radiated  from  combustion  noise  alone  through  a reflection  free 
end  plane. 

The  second  computation  involves  a correction  for  a reflecting 
termination  and  is  an  Impedance  correction.  Using  estimated,  measured 
or  calculated  isentropic  impedance  coefficients 
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is  obtained  by  setting  x i,  and  a 
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— in  the  H.  of 
Y Ap 


Eq.  (12),  H.  is  obtaind  by  setting  x * ^ and  using  the  actual  a in 

AU  Ui) 

Eqs.  (12).  nils  is  merely  an  average  correction  factor  to  account  for 

spectral  differences  in  the  acoustics  of  different  rigs.  It  is  the 

source  of  the  7,9  dB  correction  applied  in  the  previous  section.  Ihe 

product  F 11  the  amount  of  pure  combustion  noise  which  would 

eq  ter 
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be  radiated  in  the  rig  in  the  absence  of  vorticity  and  entropy  noise. 

The  third  term  is  a correction,  if  known,  for  the  spectral 
shape  of  the  combustion  noise  source,  A . It  could  be  measured  by 

U) 

the  two  microphone  technique  using  Eq.  (15).  If  the  actual  A shape 

(U 

is  known 
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It  would  actually  be  more  accurate  to  do  the  averaging  on  the 
2 2 

l**Aol  product  to  get  a product  average.  The  term 

has  actually  been  assumed  unity  in  all  comparisons  of  this  report. 

nie  fourth  and  fifth  terms  account  for  additional  noise  due  to 

entropy  and  vorticity  effects.  For  the  diffuser  runs  of  this  report 

and  considering  the  full  frequency  range  0-2500  Hz  Tl^=  1.1 , 1)^=  0.1. 

For  the  nozzle  only  run,  Tl^*  0.3,  11^=  0.2.  In  any  given  case  these 

would  be  difficult  to  estimate  but  would  really  have  to  be  measured. 

One  can  tell,  however,  whether  they  are  zero  because  there  are  some 

configurations  such  as  those  of  Refs.  (1)  and  (4)  where  they  were  zero. 

As  one  moves  to  the  right  in  Eq.  (26)  it  becomes  more  and  more 

difficult  to  estimate  the  values  of  the  terms.  71^  ^ accounts  for 

'int 

additional  noise  if  the  sources  are  correlated.  In  this  report  it  is 
conjectured  that  for  the  diffuser  runs  ea  2,  but  for  any  given 

case,  depending  on  the  type  of  correlation,  any  number  for  a three 
source  system  could  emerge  between  the  limits  0 and  4.  It  is  perhaps 
safest  to  neglect  this  effect  and  set  1. 
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'll.  accounts  for  all  effects  downstream  of  the  combustor 
uown 

termination  plane  or  the  choking  plane  In  the  case  of  a choked  nozzle 
termination.  In  the  diffuser  runs  here  w 3 If  sa  2,  as  con- 

jectured. In  the  nozzle  only  run  Tl.  __  w 1,  It  is  here  where  the 
effects  of  an  engine  installed  configuration  would  enter.  That  is, 
this  factor  would  account  for  the  acoustical  behavior  and  additional 
coherent  generation  of  noise  in  the  turbine  and  nozzle  assemblies. 

Such  considerations  are  beyond  the  scope  of  this  work. 

In  Eq,  (26)  contains  the  basic  combustor  variables  scaling 
laws.  Since  all  sources  scale  in  roughly  the  same  manner,  as  shown 

in  this  report,  P should  contain  proper  combustor  variables  scaling 
eq 

behavior.  It  is  essential  that  the  correction  be  made, and  after 

that, if  11^(1  + \ + \)  \nt  ^own  equal  to  unity,  a lower  limit 

for  P is  probably  the  result,  unless  there  is  heavy  attenuation  down- 
stream before  the  power  measurement  is  made. 

Moving  specifically  to  the  possible  corrections  to  the  prediction 
formulas  of  Refs.  (1)  and  (4),  it  is  recognized  that  there  are  two  types 
of  prediction  schemes  - one  for  combustor  rigs  and  one  for  engines.  It 
is  recommended  that  for  combustor  rigs,  where  vorticity  and  entropy 
noise  are  absent,  the  correlation  of  Eq.  (21)  with  the  numbers  of  Table 
1 be  used  for  the  sound  power  output.  This  can  be  corrected  to  the  actual 
acoustic  behavior  of  the  combustor  by  using  an  estimated  from  Figure 

1.  It  is  therefore  recommended  that  the  Pratt  and  VOiitney  and  General 
Electric  combustor  rig  formulas  be  abandoned  in  favor  of  the  one  gen- 
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Turning  to  engine  correlations  the  problem  Is  much  more  diffi- 
cult. The  Pratt  and  Whitney  method  consists  of  taking  the  combustor 
rig  formula  and  correcting  It  by  a turbine  transmission  loss  and 
nozzle  loss.  The  results  of  this  program  show,  however,  that  further 
corrections  should  be  made  for  a)  vortlclty  and  entropy  noise  b)  addi- 
tional noise  generation  by  vortlclty  and  entropy  In  the  turbine  and 
c)  combustor  acoustics.  The  first  two  are  noise  gains  and  the  third 
Is  a loss  so  that  they  are  compensating  corrections.  As  a first  ap- 
proximation, therefore,  they  may  be  expected  to  cancel  and  no  change 
Is  recommended  to  the  Pratt  & Whitney  engine  formula. 

The  General  Electric  engine  prediction  scheme  Is  not  tied  to  a 
theoretical  mechanism  as  strongly  as  Is  the  Pratt  & Whitney  scheme. 
Moreover,  the  engine  correlation  of  each  manufacturer  does  not  hold 
for  the  sound  power  results  of  each  other's  engines.  No  light  has 
been  shed  on  this  discrepancy  from  this  program.  Consequently,  for 
engine  core  noise  correlations  It  appears  that  both  correlations  must 
be  accepted  as  valid  within  the  bounds  that  the  correlations  are  only 
used  for  engines  within  the  Individual  manufactxirer's  group.  While 
this  Is  unsatisfactory,  the  current  program  has  at  least  been  success- 
ful for  giving  some  guidance  In  combustor  rig  correlations.  A future 
program  Is  needed  to  clarify  the  engine  correlation  problem. 
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IX.  Conclusions 

It  Is  imperative  in  scanning  the  conclusions  below  to  recognize 

that  some  of  them  may  be  apparatus  dependent. 

1.  A correlation,  with  a rational  theoretical  basis,  has  been 
developed  for  combustors  with  impedance  matching  termiations 
or  exhausting  directly  to  the  atmosphere.  In  such  cases  only 
direct  combustion  noise  is  present. 

2.  Due  to  the  impedance  properties  of  a choked  nozzle  there  is  a 
shift  in  radiated  combustion  noise  to  higher  frequencies  than 
would  be  expected  from  prior  data. 

3.  Contributions  to  the  radiated  noise  for  choked  nozzle  termina- 
tion have  been  shown  from  direct  combustion  and  indirect  vor- 
ticity  and  entropy  sources.  An  unambiguous  split  into  causal 
relations  between  the  noise  and  the  sources  cannot  be  made, 
however,  because  of  relatively  high  coherence  between  the 
sources. 

4.  A discrepancy  between  theoretical  and  experimental  scaling  of 
the  source  strengths  was  seen  which  serves  to  explain  the  ob- 
served sound  power  scaling  law. 

5.  When  the  choked  nozzle  is  further  terminated  by  a diffuser 
additional  coherent  noise  is  generated,  although  the  source 
is  uncertain. 


f 


Entropy  noise  is  only  important  below  about  1000  Hz,  because 
of  the  impedance  behavior  of  a choked  nozzle.  Vorticity  and 
direct  combustion  noise  are  important  over  the  entire  frequency 
spectrum  and  are  of  comparable  magnitudes. 

Measurements  of  two  acoustic  pressures  downstream  of  the  pri- 
mary combustion  zone  are  sufficient  to  isolate  the  direct  com- 
bustion noise  source  strength. 

A correlation  formula  procedure  has  been  suggested  for  combus- 
tor rig  noise  prediction.  However,  for  strict  application  it 
would  require  at  least  as  detailed  measurements  as  outlined  in 
this  report. 
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Appendix  A I 

Entropy  Tube  Results 

Introduction 

The  purpose  of  using  the  entropy  tube  facility,  shown  in  Fig,  1,  was 
to  check  the  values  of  the  admittance  and  transmittance  coefficients,  ob- 
tained from  the  hot  flow  experiments.  The  admittance  coefficients  for  a 
nozzle  are  defined  by  an  admittance  relation,  at  the  nozzle  entrance,  as 
given  by 

V+ap+ba=0  (Al) 

s s s s s 

where  V , p and  are,  respectively,  the  appropriate  transforms  of  the 
8 8 8 

plane  wave  velocity,  pressure  and  entropy  perturbations,  a and  b are, 

S 8 

respectively,  the  isentropic  and  non-isentropic  admittance  coefficients. 

Equation  (Al)  expresses  a boundary  condition  at  the  nozzle  entrance  and 

(Al) 

is  obtained  by  studying  the  unsteady  behaviour  of  the  nozzle  , The 
transmittance  coefficient  for  any  nozzle  is  defined  as  the  ratio  of  the 
appropriate  transform  of  the  pressure  pulse,  transmitted  out  of  the  nozzle, 
and  the  incident  pulse  at  the  nozzle  entrance. 

The  entropy  tube  facility  was  designed  with  the  objective  of  obtaining 
a)  a transient  isentropic  pressure  pulse  and  b)  a transient  entropy  pulse. 

By  appropriate  internal  and  external  measurements  of  pressure  fluctuations, 

:! 

and  by  using  a set  of  simplified  linear  acoustic  equations  vdiich  simulate 
the  experiments,  the  admittance  and  transmittance  coefficients  for  a given 
nozzle  are  evaluated.  These  coefficients  are  compared  with  the  correspond- 
ing values  from  the  hot  flow  experiments,  described  in  section  6 of  this  re- 
port. The  admittance  coefficients  are  also  compared  with  available  approxi- 
mate theoretical  values.  This  Appendix  describes  the  experimental  facility, 
the  theoretical  formulation,  the  data  reduction  and  analysis,  and  the  final  re- 
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suits  In  terns  of  admittance  and  transmittance  values. 

Experimental  Facility 

The  objective  was  to  obtain,  experimentally,  transient  isentropic  and 
entropy  pulses  and  record  the  nozzle  response  to  these  pulses.  With  a 
transient  pulse,  one  experimental  run  is  expected  to  yield,  through  a 
Fourier  decomposition,  all  the  desired  frequency  Information,  Of  course, 
the  highest  acceptable  frequency  will  depend,  as  explained  later  in  this 
section,  on  the  rise  time  of  the  pulse, 

A schematic  of  the  entropy  tube  facility  is  shown  in  Fig.  1,  it  consists 
of  a constant  area  two  inch  diameter  duct,  having  a uniform  flow  of  air. 

Any  test  nozzle,  whose  response  is  to  be  evaluated,  can  be  connected  to 
the  downstream  end  of  this  duct.  Airflow,  obtained  from  a 150  psia  reser- 
voir, first  passes  throu^  a pressure  regulator,  which  filters  out  the 
upstream  pressure  disturbances.  Then  the  flow  goes  through  a porous  plate 
and  a rubber  hose,  both  of  which  help  to  reduce  the  flow  noise,  A static 
pressure  measuring  device  and  a flow  meter  are  provided  to  monitor  the 
steady  state  flow  variables. 

Two  independent  techniques  were  used  for  generating  the  desired  pulses. 
The  first  one,  shown  in  Fig,  1,  consisted  of  sending  a shot  of  freon,  a high 
molecular  weight  gas,  into  a uniform  airflow,  A solenoid  valve  was  used  for 
creating  the  sharp  transient  pulses.  When  switched  on,  its  operation  is 
analogous  to  the  action  of  puncturing  a diaphragm  in  a shock  tube  test.  The 
process  initiates  a series  of  left  running  expansion  pulses  and  a right  run- 
ning compression  pulse,  followed  by  the  freon  flow.  The  right  running  com- 
pression pulse,  termed  as  the  isentropic  (or  initial)  pressure  pulse,  was 
used  to  obtain  the  isentropic  admittance  and  transmittance  coefficients  of 
the  nozzle.  The  freon  front  travels  through  the  injector,  into  the  constant 
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area  duct,  vdiere  It  mixes  with  the  air  to  generate  a sharp  entropy  pulse. 

Ihls  entropy  pulse  Interacts  with  the  nozzle  to  generate  a sharp  pressure 
(A2) 

pulse'  , This  pressure  pulse,  termed  as  the  entropy  pressure  pulse,  was  used 
to  evaluate  the  non-lsentropic  admittance  and  transmittance  coefficients. 

The  typical  Isentroplc  and  entropy  pressure  pulses,  as  monitored  by  two  ml- 
cro^ones  in  the  duct,  are  shown  In  Fig.  2. 

Ihe  second  technique,  shown  In  Fig.  3,  consisted  of  creating  an  abrupt 
spark  by  discharging  a high  charge  across  two  electrodes.  The  transient 
spark  creates  an  initial  pressure  pulse  (due  to  the  energy  added)  and  an 
entropy  pulse  (or  the  hot  spot).  These  pulses  were  used,  as  in  the  freon 
injection  case,  for  estimating  the  isentropic  and  non-isentropic  acoustical 
properties  of  a given  nozzle.  The  typical  Isentropic  and  entropy  pressure 
pulses  for  the  spark  case,  as  monitored  by  two  microphones  in  the  duct,  are 


shown  in  Fig,  4. 

A comparison  of  Figs.  2 and  4 shows  the  essential  difference  between 
the  pulses,  generated  by  the  above  two  techniques.  The  Initial  pressure 
pulses  look  very  similar  except  that  the  pulse  generated  by  the  spark  has 
a much  higher  amplitude  as  compared  to  the  one  generated  by  freon  injection. 

Also  the  spark  creates  a much  sharper  initial  pulse,  i.e.  the  rise  time  of 

the  pulse,  is  much  smaller.  Some  explanation  about  the  entropy  pressure  ' > * 


pulses,  generated  by  these  two  techniques,  seems  appropriate  here.  The 
freon  Injection  process  creates  a high  molecular  wel^t  or  a positive  den- 
sity pulse,  which,  encountering  the  nozzle,  generates  a compression  pulse. 
This  is  shown  as  pulse  (2)  in  Fig,  2,  On  the  other  hand,  the  spark  discharge 
is  equivalent  to  a negative  density  pulse  and  thus  results  in  an  expansion 
pulse.  (These  conclusions  are  theoretically  substantiated  in  Ref.  (A2)). 
However,  the  theoretical  formulation,  given  in  a following  section,  has 


TIME  (MILLISECONDS) 

FIGURE  2,  TYPICAL  SIGNALS  AS  RECORDED  BY  MICROPHONES  Pi  AND  P2  IN  THE 
ENTROPY  TUBE.  FREON  INJECTION  CASE. 


FIGURE  4,  TYPICAL  SIGNALS  AS  RECORDED  BY  MICROPHONES  P|  AND  Po  IN  THE  ENTROPY  TUBE. 
SPARK  CASE. 


been  made  independent  of  the  nature  of  the  entropy  pressure  pulse.  Thus, 
both  techniques  should  be  Identical  as  far  as  the  evaluation  of  nozzle 
properties  is  concerned. 

For  monitoring  the  pressure  pulses  inside  the  tube,  microphones  with 
an  infinite  tube  arrangement,  as  shown  in  Fig.  1,  were  used.  With  the  in- 
finite tube,  the  microphones  are  not  in  direct  contact  with  the  flow  field 

and  thus  the  non-propagational  hydrodynamic  noise  contamination  is  reduced 

(A3) 

a substantial  amount^  ' » The  conical  nozzle  used  for  the  present  program 
is  shown  in  Fig.  5.  Ideally  this  nozzle  should  be  designed  for  the  flow 
conditions  used  in  the  how  flow  experiments,  where  the  upstream  Mach  num- 
ber was  about  .02.  But  the  lowest  possible  Mach  number  in  the  entropy  tube 
facility  was  about  .06.  Hence  the  nozzle  for  the  present  program  was  designed 
for  an  entrance  Mach  number  of  .065  and  choked  conditions  at  the  throat.  The 
length  of  the  nozzle  was  selected  to  give  the  same  Mach  rnanber  gradient  as 

used  in  the  hot  flow  studies,  since  it  is  this  gradient  that  is  responsible 

(Al) 

for  noise  generation  by  incident  entropy  waves  . 

It  was  observed  that  for  choked  exhaust  conditions,  the  Jet  noise  was 
blanketing  the  pressure  pulses  transmitted  out  of  the  nozzle.  To  alleviate 
this  difficulty,  a diffuser  was  installed  at  the  downstream  end  of  the  nozzle. 
The  diffuser  was  designed  to  reduce  the  exhaust  velocity  to  an  insignificant 
value.  The  5°  diffuser  angle  was  selected  to  avoid  stall  possibilities.  In 
fact,  the  nozzle  diffuser  assembly  served  two  important  objectives.  It  a) 
reduced  the  Jet  noise  contamination  significantly  and  b)  preserved  choked 
nozzle  boundary  condition:  This  condition  is  important  if  the  experimental 
results  are  to  be  compared  with  the  approximate  theoretical  results,  avail- 
able only  for  choked  conditions.  Internal  pressure  measurements  were  used 
for  evaluating  the  isentropic  and  non-isentropic  admittance  values  of  the 


nozzle.  The  external  pressure  measurements  were  used  for  obtaining  the 


nozzle  transmittance  coefficients. 


Theory 


The  development  of  the  governing  eqtiations,  which  simulate  the  experi- 
mental conditions,  is  deleted  here.  Only  the  final  resulting  equations, 
which  make  use  of  the  experimental  data,  are  presented  (For  details,  see 
Ref.  (A2)).The  isentropic  admittance  coefficient  a and  the  reflection 

coefficient  R of  the  terminating  nozzle  are  given  by 
s 


s 


1-2 


gi(mj^-hii2)(Z2-Zi)  _ ^ 


ei°*l(Z2~Zl)  Ps^/Psa] 
|Psi/pS2)  - e“^2(Z2-Zl) 


,-i(ii^-hn2)Z2 


where  is  the  ratio  of  the  specific  heats  of  air,  Mq  is  the  flow  Mach 
number  at  the  nozzle  entrance,  Z^  and  Z2  are  the  dimensionless  (diameter 
of  Che  tube  used  for  non-dimensionallzaClon)  locations  of  Che  microphones 
1 and  2 from  the  source  position,  mj  and  m2,  which  represent  the  wave 
nimbers  for  the  right  running  and  left  running  waves,  are  given  by 

Pl“o 
“1  “ rTM; 

(A4) 

^1^ 


where  is  the  frequency,  non-dimensionalized  by  using  Che  duct  diameter 
and  the  steady  state  velocity  in  Che  duct.  Finally,  in  Eqs.  (A2)  and  (A3), 
Pgj^  and  Pg^  the  Laplace  transforms  of  the  isentropic  pressure  pulse, 
as  measured  by  the  microphones  1 and  2.  It  is  worth  mentioning  here  that 
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the  Laplace  transform  (LT)  vas  used  in  the  present  analysis.  To  get  the 
frequency  Infoxnation  for  a non-perlodlc  signal,  as  is  the  present  case, 
the  normal  procedure  is  to  make  use  of  the  Fourier  transform  (FT)  technique. 
But  for  the  signals  generated  in  the  present  experiments,  the  FT  does  not 
exist;  only  the  LT  exists.  Therefore,  a method  was  devised  to  obtain  the 
LT  on  the  Fourier  analyser.  The  details  about  the  non-existence  of  FT  and 
the  method  of  obtaining  LT  are  given  in  Ref.  (A2). 

According  to  Eqs.  (A2)  and  (A3),  the  isentropic  admittance  coefficient 
and  the  reflection  coefficient  may  be  evaluated  by  obtaining  the  Laplace 
transforms  of  the  initial  pulse,  as  recorded  by  two  microphones. 

The  entropy  pressure  pulse,  generated  as  a result  of  a cold  (or  hot) 
spot  interacting  with  the  pressure  gradient  in  the  nozzle,  was  used  to 
estimate  the  non-iaentropic  admittance  coefficient.  The  LT  of  this  pulse 
is  given  by 


b Y M 
s a o 

?8a  “ 1 -vAag 


(A5) 


where  bg  is  the  non-isentropic  admittance  coefficient.  Qg  represents  the 
strength  of  cold  (or  hot)  spot.  The  attempt  to  measure  Og  directly  had 
to  be  abondoned  because  of  several  practical  problems.  Instead,  the  eval- 
uation of  be  was  carried  out  in  an  indirect  way.  It  was  assumed  that 
could  be  maintained  identical  for  runs  with  two  different  nozzles.  This 
assumption  was  completely  valid  for  the  spark  case,  as  all  the  source 
variables  could  be  maintained  constant.  For  the  freon  case,  the  freon  in- 
jection process,  for  two  runs  with  different  nozzles,  was  compared  by 
monitoring  the  pressure  pulses  in  the  freon  line.  Figure  6 gives  this 
comparison  for  the  experimental  conditions  used  in  Ref.  (A2),  whereas 


I; 

i* 
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COMPARISON  OF  FREON  LINE  SIGNALS  FOR  TWO  RUNS  WITH  DIFFERENT  NOZZLES 
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Fig,  7 gives  the  same  for  the  present  experiments.  According  to  Fig.  6,  the 
above  assumption  Is  perfectly  valid  for  the  experiments  reported  in  Ref,  (A2), 
However  for  the  present  experiments,  as  Fig,  7 indicates,  the  validity  of 
the  above  assumption  could  not  be  completely  established.  The  reason  is  dis- 
cussed below. 

In  order  to  obtain  in  the  present  experiments,  one  run  was  made  with 
an  orifice  plate,  which  simulates  the  quasi-steady  or  zero  frequency  case, 
and  the  other  one  with  the  conical  test  nozzle.  The  exhaust  velocity 
with  the  orifice  plate  was  about  1000  f/s  as  compared  to  a value  of 
70  f/s  for  the  conical  nozzle,  as  the  conical  nozzle  had  a diffuser  at  its 
downstream  end.  This  amounted  to  a much  higher  static  pressure  in  the  en- 
tropy tube  for  the  orifice  plate  as  compared  to  that  for  the  conical  nozzle, 
thereby  causing  unequal  pressure  differences  across  the  solenoid  for  the  two 
runs.  Some  corrections,  as  explained  in  a later  section,  were  incorporated 
to  offset  these  different  experimental  conditions.  Then  with  the  above 
assumption,  Eq.  (A5)  leads  to  the  following  expression  for  the  non- 
Isentroplc  admittance  coefficient  (bs^)  for  the  test  nozzle. 

•’so  Usjjjo  “o  ®So  (A6) 


where  a and  bg  are  the  theoretical  admittance  coefficients  at  zero  fre- 
o o 


quency  and  are  given  by 


(A7) 
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and 


TYPICAL  FREON  LINE  MICROPHONES'  SIGNALS 


FIGURE  7.  COMPARISON  OF  FREON  LINE  SIGNALS  FOR  RUNS  WIIHIIHE  CONICAL  TEST  NOZZLE 
AND  AN  ORIFICE  PIATE. 


Thus  the  non-lsentroplc  admittance  coefficient  may  be  evaluated  by  measuring 
the  Laplace  transform  of  the  entropy  pressure  pulses  for  the  orifice  plate 
and  the  test-nozzle. 

The  transmittance  coefficient  or  the  transfer  function  relates  the  out- 
put and  input  signals  for  a given  nozzle.  It  may  be  defined  as  the  ratio 
of  the  appropriate  transforms  of  output  and  input  signals.  For  the  isen- 
tropic  transmittance  coefficient  the  output  is  the  isentropic  pressure  pulse, 
transmitted  out  of  nozzle,  whereas  the  input  is  the  incident  pressure  pulse. 
Mathematically,  it  may  be  written  as 


where  p is  the  transform  of  the  transmitted  Isentropic  pulse  and  Pg  is 

the  transform  of  the  incident  pressure  pulse.  Now  the  pulse,  monitored  by 

the  microphone  at  the  nozzle  entrance,  contains  both  the  Incident  pressure 

pulse  (p_  ) as  well  as  the  reflected  pulse  (p„  ),  i.e. 

®1  r 


Psi  + Ps, 


where  Pg  and  pg  are  related  through  Rg  by  the  equation 

Pb  ” K Pb  . 


(AlO) 


Use  of  Eqs.  (A9)  and  (AlO)  in  (A8)  gives 

Pti 

T„  “ — ^ ( 1 + R ) 

«i  Ps2  ® 


(All) 


t 


Thus  the  isentropic  transmittance  coefficient  may  be  evaluated  by  the 
appropriate  transforms  of  the  initial  pulse,  as  monitored  by  the  various 
microphones,  in  combination  with  Rg,  which  is  given  by  Eq.  (A3). 
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For  the  non-isentropic  transmittance  coefficient  (Tg  ),  the  output  is 

the  transmitted  non-isentropic  pressure  pulse,  \»hereas  the  input  is  given 

by  o_»  the  strength  of  the  cold  (or  hot)  spot.  With  the  aid  of  Eq.  (A5), 
s 

Tq  can  be  written  as 


Pi-  bq  V M 

p„  1 - Y a„ 

^ o St 


(A12) 


where  Pj^  is  the  transform  of  the  transmitted  entropy  pressure  pulse.  The 
evaluation  of  b„  and  a^  for  the  test  nozzle  was  discussed  above.  Then 

St  St 

using  the  computed  ag^  and  bg^  values  and  the  transforms  of  the  entropy 
pressure  pulse,  the  non-isentropic  transmittance  coefficient  may  be  eval- 
uated. 

Data  Reduction  and  Analysis 

The  sequence  of  the  experimental  and  theoretical  steps  leading  to 
the  evaluation  of  the  acoustical  properties  of  a given  nozzle  is  given  in 
Figs,  8 and  9.  Figure  8 illustrates  the  experimental  procedure  leading  to 
the  evaluation  of  ag,  Rg,  b^,  Tg^  and  Tg^,  As  shown  in  Fig,  8,  the  required 
transient  pulses  are  generated  either  by  freon  injection  or  by  spark  dis- 
charge. These  pulses,  inside  the  tube  and  outside  in  an  anechoic  chamber, 
are  monitored  by  various  microphones  and  recorded  on  magnetic  tape.  The 
recorded  data  is  digitised  and  properly  averaged  in  the  Fourier  analyser 
(see  Ref.  (A2)  for  details).  After  separating  the  isentropic  and  the  non- 
isentropic  pulses,  their  Laplace  transforms  (LT)  are  obtained.  These  LT 
values  are  used  in  Eqs,  (A2),  (A3)  and  (A6)  to  yield  the  nozzle  admittance 
coefficients.  Similarly,  the  LT  of  the  isentropic  and  non-isentropic  trans- 
mitted pulses,  in  combination  with  the  derived  nozzle  admittance  values,  is 
used  in  Eqs,  (All)  and  (A12)  to  obtain  the  transmittance  coefficients. 
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Figure  9 lists,  briefly,  the  steps  involved  in  the  theoretical  formula- 
tion, which  uses  the  experimental  results  to  yield  the  desired  nozzle  proper- 
ties. The  data  acquisition,  reduction  and  analysis  procedure  is  illustrated 
in  Fig,  10. 


Results  and  Discussion 

The  results  obtained  by  using  the  experimental  data  and  the  theoretical 
relations  are  shown  in  Figs,  11-23,  Wherever  possible,  the  corresponding 
theoretical  values,  obtained  independently,  are  also  plotted.  For  these 
plots,  a dimensionless  frequency  (g)  has  been  used.  It  is  defined  by  the 
relation 


du/dx 


(A13) 


where  uj  is  the  frequency  in  radians  per  second  and  ^ is  the  steady  state 
velocity  gradient  for  a given  nozzle.  Since  ~ is,  in  general,  not  constant 
in  the  nozzle,  only  an  overall  or  an  effective  ^ can  be  used  in  Eq,  (A13). 
The  importance  of  using  P becomes  obvious  if  it  is  mentioned  that  the  quasi- 
steady or  zero  frequency  case  corresponds  to  ^"*0,  rather  than  Essen- 

tially, the  velocity  gradient,  an  Important  parameter  in  the  noise  response 
studies,  is  thus  incorporated  in  the  above  definition. 

In  the  theoretical  evaluation  of  nozzle  admittance  coefficients  (see 
Refs.  (A2)),  du/dx  is  assumed  to  be  a constant.  A nozzle,  with  a constant 
^ ,is  termed  here  as  a linear  nozzle.  This  simplifies  the  theoretical  eval- 
uation '•.onsiderably.  For  a linear  nozzle,  Eq,  (A13)  reduces  to 


a (A14) 

where  a is  the  speed  of  sound  at  the  nozzle  throat,  x is  the  distance  be- 
tween the  throat  and  the  location  of  zero  u.  This  location  is  easily  calcul- 
able for  a linear  nozzle. 
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KGURE  10,  SCHEMATIC  OF  THE  DATA  ACQUISITION  AND 
REDUCTION  PROCEDURE. 
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In  the  entropy  tube  experiments,  it  was  required  to  simulate  the  con- 
ditions of  the  hot  flow  experiments,  where  a conical  nozzle  was  employed. 

* 

For  the  conical  nozzle,  it  was  assumed,  as  an  approximation,  that  x rep- 
resents the  nozzle  length,  lliis,  in  fact,  defines  an  effective  du/dx  in 
Eq.  (A13).  Ihen,  for  the  hot  flow  experiments,  where  the  results  are  ob- 
tained for  a frequency  range  of  0-2000  Hz,  the  p values  range  from  0 to  1.8. 
Before  discussing  any  results, it  should  be  pointed  out  again  that  for  the 
entropy  tube  experiments,  a conical  nozzle  was  used,  whereas  the  theoretical 
results  were  obtained  for  a linear  nozzle.  Whenever  the  experimental  results 
are  compared  with  the  corresponding  theoretical  values,  this  difference  in 
the  experimental  and  theoretical  nozzlea  must  be  borne  in  mind. 


Results  From  Freon  Injection  Experiments 

As  mentioned  earlier,  the  evaluation  of  Psj^/Ps2  determine  a^  and 

Rg,  Hence,  P8i/Ps2  * experimentally  as  a function  of  (Hz)  and 

P,  is  plotted  in  Fig.  11,  Also  shown  in  this  figure,  is  the  theoretically 

estimated  Pg2^/Pg2»  (the  procedure  to  obtain  theoretical  Pgj/Ps2  described 

in  Ref,  (A2)).  According  to  Fig,  11,  the  experimental  and  the  theoretical 

Po, /Po«  agree  reasonably  well  up  to  a frequency  of  225  Hz,  This  corresponds 
®1  “2 

to  a p value  of  about  .50.  This  frequency  limit  was  imposed  by  the  nature 
of  the  Initial  pressure  pulse,  generated  by  the  freon  injection  apparatus. 
For  frequencies  higher  than  225  Hz,  the  smallest  resolvable  value  in  the 
Fourier  analyser  dominated  the  initial  pulse  components.  Both  the  experi- 
mental and  the  theoretical  ag  values  are  plotted  In  Fig.  12.  The  correspond- 
ing reflection  coefficients  are  shown  in  Fig,  13,  The  amplitude  of  Sg,  ob- 
tained experimentally,  shows  an  acceptable  agreement  with  the  theory  for 
frequencies  between  100  and  250  Hz.  The  reason  for  an  unacceptable  variance 
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between  100  and  250  Hz,  the  quantitative  agreement  between  experiment  and 
theory  is  poor;  the  reason  being  the  — factor  in  Eq.  (A2),  For  the 
present  experiments,  this  factor  has  a value  of  about  116,  which  means  that 


a small  error  in  pg^^/pg^svaluation  gets  multiplied  by  a factor  of  about  116, 
thereby  causing  a poor  quantitative  agreement.  However,  the  a^  phase  is  not 
affected  by  this  factor,  and  thus  shows  a better  agreement  with  theory.  The 
conclusions  are  further  supported  by  the  plots  in  Fig,  13,  As  the  error 
magnification  reasoning  does  not  apply  to  the  Rg  evaluation.  Fig,  13  shows 
a much  better  accord  between  experiment  and  theory.  As  mentioned  in  the  be- 
ginning of  this  section,  another  possible  cause  for  the  disagreement,  be- 
tween the  experiment  and  theory,  could  be  the  fact  that  the  theoretical  and 
experimental  nozzles  were  not  matched  exactly.  It  is  not  possible,  with  the 
theory  used,  to  predict  the  error  caused  by  this  mismatch.  However,  in  Ref. 
(A2),  where  the  nozzles  were  matched,  better  results  were  obtained. 

To  estimate  the  non-isentropic  admittance  coefficient,  freon  injection 

runs  were  made  with  the  orifice  plate  and  the  test  nozzle.  The  corresponding 

(Pa  ) and  (p-  )_,  thus  evaluated,  were  used  in  Eq,  (A6)  to  obtain  bg  /bg  , 
O X u o 

At  low  frequencies,  and  Sg^  (or  ago)  are  both  of  the  order  of  ,1.  Thus  the 
term  Ya*^®8t  ^ ®So^»  small  compared  to  unity,  was  neglected 

in  Eq,  (A6).  As  mentioned  in  a previous  section,  the  underlying  assumption 
(that  a is  same  for  runs  with  different  nozzles)  could  not  be  completely 
satisfied  in  the  present  experiments.  To  offset  the  effect  of  unequal  pres- 
sure differences  across  the  solenoid  for  the  orifice  plate  and  the  conical 


I 


118 


nozzle,  some  corrections  were  Included.  These  corrections  were  based  on 
the  shock  tube  relations,  as  given  In  Ref.  (A4).  The  final  corrected  ex- 
perimental bg^/bs^  as  well  as  the  corresponding  theoretical  values  are 
shown  In  Fig.  14.  According  to  this  figure,  the  bg^  values  for  the  conical 
nozzle  are  reasonably  acceptable  up  to  a frequency  of  175  Hz  (8^.4).  Above 
this  frequency,  the  smallest  resolvable  value  in  the  Fourier  anayser  again 
dominated  the  entropy  pressure  pulse  components. 

Next,  by  using  the  experimental  Ptj/Ps2  theoretical  Rg  values 

in  Eq.  (All),  the  Isentroplc  transmittance  coefficient  was  obtained.  As 
shown  in  Fig.  15,  the  amplitude  of  Tg^  increases  with  frequency, 
except  between  150-250  Hz,  where  it  seems  to  go  through  a minimum.  This 
can  be  explained  by  a closer  look  at  Eq.  (All)  and  Fig.  (13).  In  the  range 
150-250  Hz,  Rg  has  a phase  of  about  180°,  which  means  that  the  term  (1+Rg) 
in  Eq.  (All)  is  close  to  zero.  Consequently,  the  T^^  magnitude  is  small. 
Acttially,  both  (1  + Rs)  and  pg2  should  have  been  close  to  zero,  thereby 
giving  a 0/0  form  to  the  right  hand  side  of  Eq.  (A. 11).  But  the  error  in 


the  experimental  evaluation  of  pg^  (though  Pg2  does  have  a small  value,  as 
indicated  by  Fig.  16)  makes  the  value  of  ^*^i  to  be  finite  and  thus  causes 


T to  be  small.  Moreover  it  is  important  to  point  out  that  the  use  of  the 
®i 

experimental  Rg  Instead  of  the  theoretical  one  seems  to  yield,  as  shown  in 


Fig,  17,  much  more  acceptable  Tg^  values  in  the  troublesome  frequency  range. 
There  is  nothing  strange  about  this  behavior;  the  errors  in  the  experimental 


I 
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evaluation  of  (1  + Rg)  and  pg2  possibly  cause  the  amplitudes  of  these  two 
terms  to  be  non-zero  and  thus  improve  the  Tg^  evaluation  in  that  frequency 
range. 

The  results  for  the  non-isentropic  transmittance  coefficient  Tg^,  for  the 
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FIGURE  14,  NON-ISENTRDPIC  ADMITTANCE  COEFFICIENT  (bg  /bg  ) FOR  THE  FREON 
INJECTION  CASE.  ^ ° 


■ISENTROPIC  TRANSIITTANCE  COEFFICIENT  (To  ) FOR  FREON  INJECTION  CASE 
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mental  p»,  /pa  and  the  theoretical  a.^  and  b_^  values  (the  reason  for  not 
Rj  ®t  “t 

using  the  experimental  ag^  and  bg^  is  discussed  in  Kef,  (A2))  in  Eq.  (A12), 
According  to  Fig,  18,  the  results  seem  to  be  acceptable  up  to  a frequency 
of  175  Hz  M .4),  which  is  the  upper  limiting  value  for  the  entropy  pres- 
sure pulse,  generated  by  the  freon  injection  apparatus. 

Spark  Discharge  Data; 

The  freon  injection  data  had  an  upper  frequency  limit  of  225  Hz  for  the 
isentropic  pulse  and  175  Hz  for  the  non-lsentropic  pulse.  This  frequency 
limit  is  apparatus  dependent,  i.e,  with  a reduction  in  the  rise  time  of  the 
pulse,  it  is  possible  to  raise  the  frequency  limit.  The  spark  discharge  pro- 
vides a good  illustration.  The  isentropic  pressure  pulse  generated  by  the 
spark  discharge  has  a shorter  rise  time  as  compared  to  the  one  obtained 
by  freon  injection.  Consequently,  the  upper  frequency  limit  for  the  spark 

case  is  about  350  Hz  (p  w ,75),  Therefore,  as  shown  in  Figs,  19-21,  theory 

p 

and  experiment,  for  ^1  , ag  and  evaluations,  agree  reasonably  well  up 

Ps2 

to  350  Hz,  Again,  as  in  the  freon  case,  at  frequencies  close  to  zero,  the 
experimental  a^  values  are  not  acceptable  because  of  the  ^ form  of  the  right 
hand  side  of  Eq,  (A2),  Also,  between  50-350  Hz,  the  experimental  Sg  is  con- 
sistently higher  than  its  theoretical  counterpart,  because  the  term  , 

in  Eq,  (A2),  has  a large  value.  This  has  also  been  pointed  out,  when  the  freon 
results  were  discussed.  Furthermore,  in  Ref,  (A2),  where  experimental  and 
theoretical  conditions  were  matched  exactly,  the  spark  discharge  method 
yielded  better  results. 

Following  the  same  procedure  as  used  in  the  freon  injection  case,  the 
isentropic  transmittance  coefficient  Tg^  was  computed.  Figure  (22)  shows 
the  magnitude  and  phase  of  Tgj^,  obtained  by  using  the  theoretical  Rg  values. 


r- 


FIGURE  19.  RATIO  OF  THE  lAPIACE  TRANSFORMS  OF  THE  ISENTROPIC  PRESSURE  PULSE 
AS  MEASURED  BY  TWO  DUCT  MICROPHONES  FOR  SPARK  CASE. 


FIGURE  20.  ISENTROPIC  ADMITIANCE  COEFFICIENT  (a.)  FOR  THE  SPARK 


The  upper  acceptable  frequency  for  with  spark  discharge,  is  about  500  Hz 
> 1.1)  as  compared  to  350  Hz  with  freon  injection.  Again,  as  shown  in  Fig, 
(22),  the  Tg^  magnitude  becomes  small  at  frequency  values  of  100,  300  and  500 
Hz,  where  the  Rg  phase  is  close  to  180°,  The  reason  for  the  objectionably  small 
Tg^  values  at  these  frequency  points  becomes  obvious  if  one  follows  the  dis- 
cussion of  the  preceding  section.  Moreover,  as  shown  in  Fig.  (23),  an  improve- 
ment in  the  T values  at  these  three  points  seems  possible  if  the  experimental 
®i 

is  preferred  over  its  theoretical  counterpart. 

The  estimates  of  non-isentropic  admittance  and  transmittance  coefficients 

by  using  the  spark  discharge  method  are  not  included  in  this  report.  A few 

trial  runs  were  made  and  the  ^St^^So  values,  obtained  by  the  procedure  discussed 

in  a previous  section,  were  good  up  to  about  150  Hz.  This  frequency  limit  was 

no  better  than  the  one  obtained  by  the  freon  injection  process.  Moreover,  for 

the  spark  case,  some  of  the  uncertainties,  regarding  the  Identification  of 

the  entropy  pressure  pulse  and  the  repeatability  of  data,  could  not  be  resolved 

within  the  time  limit  of  this  program.  Consequently,  the  b.^  and  1.  evalua- 

“ t o 

tions  for  the  spark  discharge  case  are  omitted  from  this  report. 
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Appendix  B 

Cooled  Film  Probe  Results 
Basic  Theory  of  Cooled  Film  Measurements 

The  main  objective  of  the  cooled  film  measurements  Is  to  obtain 
the  Impedance  condition  at  the  nozzle  entrance.  This  requires  the 
knowledge  of  pressure,  temperature  and  velocity  fluctuation  details. 

The  pressure  and  temperature  fluctuations  were  measured  through  Photo- 
con  transducers  and  thermocouples,  whereas  the  velocity  fluctuations 
were  to  be  obtained  through  cooled  film  measurements. 

The  cooled  film  anemometer  Is  a device  to  make  velocity  fluctua- 
tion measurements  in  high  temperature  environments  similar  to  those 
obtainable  at  normal  temperatures  with  a hot  wire  or  hot  film  ane- 
mometers. The  normal  hot  wire  or  film  is  limited  by  the  maximum  op- 
erating temperature  of  the  sensor.  Since  the  actual  measurement  de- 
pends on  the  heat  transfer  mechanism,  a non-cooled  sensor  must  be 
operated  at  higher  than  environmental  temperatures.  The  cooled  sensor 
circumvents  this  problem  by  adding  an  Internal  heat  sink  so  that  Its 
surface  Is  actually  below  the  temperature  of  the  environment.  In  this 
manner  It  Is  possible  to  expose  sensors  to  very  high  temperature  en- 
vironments, the  only  limitation  being  the  maximum  cooling  rate  that 
can  be  applied  to  the  sensor. 

The  sensing  element  consists  of  a pyrex-u-tube,  15  mm  diameter, 
through  which  a cooling  fluid  can  be  circulated.  On  one  leg  of  this 
U-tube  a platinum  film  of  one  micron  thickness  and  1.5  mm  length  Is 
deposited  and  a gold  coating  on  the  rest  of  the  tube  serves  as  the 
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electrical  contact.  With  the  coolant  circulating  through  the  sensor, 
the  platinum  film  Is  maintained  at  constant  reslsLunce  (and  hence 
constant  temperature)  by  means  of  a bridge  circuit,  of  which  the 
sensor  forms  a part, and  a feedback  system  which  utilizes  voltage 
and  power  amplifiers  to  maintain  bridge  balance.  According  to  its 
manufacturers  (TSl),  the  frequency  response  of  this  heat  flux  system 
is  rated  at  50  kHz. 

Figure  1 shows  the  arrangement  of  the  cooled  film  sensor  and 
Che  coolant  system.  Hie  basic  heat  transfer  mechanism  involved  in 
the  cooled  film  measurements  is  illustrated  in  Fig.  2.  The  heat 
transfer  relation  can  be  written  as, 

P + Qe  - (1) 

where  E 2 

P ■ electrical  power  input  to  the  sensor  = ( ) R (2) 

E ■ Bridge  voltage 

R =■  Sensor  resistance 

Q > Heat  transferred  from  the  sensor  surface  to  Che 
cooling  fluid 

Heat  transferred  from  the  environment  to  the 
sensor  surface.  (2) 

In  Che  absence  of  any  heat  transfer  between  the  film  and  Che  environ- 
ment, the  electrical  power  supplied  to  the  film  element  is  dissipated 
by  the  coolant  circulating  through  the  sensor.  When  this  element  is 
exposed  to  environmental  conditions  such  that  a power  transfer  be- 
tween the  environment  and  the  element  occurs,  the  power  supplied  to 
the  element  by  the  heat  flux  system  is  adjusted  so  that  the  power 
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CONNECTION  TO  SENSOR 


Figure  1.  (k>oled  Film  Arrangement. 
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dissipation  to  the  coolant  is  unaltered.  In  this  way,  the  fluctua- 
tions in  the  bridge  voltage  are  indicative  of  changes  in  the  environ- 
mental conditions. 

The  cooled  film  sensor  not  only  registers  velocity  fluctuations 
but  also  is  sensitive  to  pressure  and  temperature  fluctuations.  Hence, 
the  primary  task  is  to  relate  the  cooled  film  output  to  these  fluctua- 
tions. This  Involves  an  empirical  heat  transfer  relation  by  Ahmed 
This  correlation  is  given  by 


c + d (Re^)" 


(3) 


where  istheNusselt  Number,  is  the  Reynolds  Number,  v is  the 
kinematic  viscosity,  c,  d,  m & n are  empirical  constants  and  the  sub- 
scripts e and  f refer  to  environmental  and  film  conditions.  Writing 
Nusse It  Numbers  intermsof  heat  transfer  rates  and  after  substituting 
Eqs.  (1)  and  (2)  into  Eq.  (3),  the  resulting  cooled  film  voltage  fluc- 
tuations can  be  written  as 


E' 

g 


+ B 


(4) 


where  A,  B and  C depend  on  the  operating  conditions  for  a typical 
case.  In  the  present  program,  A ■ - 0.7  B = -0.7  and  C = -4.79. 


Admittance  Calculations 

The  admittance  condition  at  the  entrance  to  the  nozzle,  after 


i 

! ‘ 


averaging  over  cross-sectional  area,  is  written  as 


T 

A = cross  sectional  area 
e 

a^>  = Isentropic  and  Entropy  admittance  coefficients 


The  cooled  film  relation  of  Eq,  (4),  after  a Fourier  transform  opera- 
tion and  area  averaging,  can  be  written  as 

^ " A J ^^^o  O) 

Substitution  of  Eqs.  (6)  and  (7)  into  Eq.  (5)  yields 

R + a p^  + b J = 0 (8) 

W O OP  ' ' 

where  a * a B - b B - A (9) 

(U  U)  U>  Y 

b - b B - C 
U)  U) 

and  J ■ T f tdA 

A ^ 

Equation(8)contalns  two  unknowns  a^  and  b^.  By  performing  cross  spectral 
operations  with  any  two  of  the  quantities  R,  p^  and  J,  two  equations 
can  be  obtained  for  two  unknowns.  The  equations  used  in  the  present 
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program  are 


S„  + aS  +bS=0 
Pn^  OUftjIb  HJfbJ 


S„  + a S.  +b  S„  = 0 
JR  Pq  U)  JJ 


(10) 


where  S refers  to  the  cross  spectral  operation  between  the  sub- 
scripted quantities.  It  is  important  to  note  that  S and  S in 
Eq,  (10),  have  been  obtained  after  converting  the  single  point  meas- 
urements to  an  area  averaged  one,  through  the  use  of  correlation 

/<> 

area  values.  Equation(10)can  be  solved  for  a and  b and  the  im- 

U)  (JU 

pedance  coefficients  a and  b are  then  obtained  through  Eqs.  (9). 

(JU  U) 


Results 

The  manufacturer's  claims  for  this  probe  were  shown  to  be  totally 
false.  From  Eq,  (4)  and  the  stated  A,  B,  and  C coefficients  it  is 
seen  that  the  probe  acts  best  as  a thermometer.  Moreover,  the  com- 
bustor has  extremely  high  temperature  fluctuations  so  that  a meas- 
urement was  taken  of  a probe  and  a thermocouple  placed  very  close  to 
each  other  in  the  exhaust  of  the  combustor.  If  the  C values  was  in- 
deed constant  and  the  thermocouple  was  compensated,  as  it  was,  the 
transfer  function  between  the  two  measurements  should  have  been -4.8. 
Shown  in  Fig.  3 is  the  transfer  function.  Clearly^ a)  even  at  low  fre- 
quency the  transfer  function,  while  negative,  is  not  large  and  b) 
serious  frequency  response  starts  at  around  80  Hz. 

Data  reduction  for  the  admittance  coefficients  was  nevertheless 


attempted  on  a run  with  the  admittance  results  of  Fig.  4.  The  b^  is 


TRANSFER  FUNCTION  (normolized) 


Figure  3.  Transfer  function  between  the  cooled  film  probe  and 
thermocouple. 


MITTANCE  CQEFFICI 


Figure  4.  Admittance  coefficients  derived  from  the  cooled  film  probe 
results  compared  with  the  theoretical  results. 
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off  by  a factor  of  at  least  10  when  compared  with  tinory  and  there 
la  erratic  behavior  of  the  Isentroplc  admittance  coefficient. 

vnille  some  work  was  done  on  the  compensation  of  the  probe,  it 
soon  become  evident  that  substantial  work  would  be  required  to  fully 
characterize  Its  behavior.  Consequently,  the  use  of  the  probe  in  this 
program  was  abandoned.  Ihe  probe  gives  a good  signal  which  should  be 
useful  in  future  work,  but  it  needs  full  characterization. 
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